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In these Canberra summer school lectures we treat a number of topical issues 
in neutrino astrophysics: the solar neutrino problem, including the physics of the 
standard solar model, helioseismology, and the possibility that the solution involves 
new particle physics; atmospheric neutrinos; Dirac and Majorana neutrino masses 
and their consequences for low-energy weak interactions; red giant evolution as a 
test of new particle astrophysics; the supernova mechanism; spin-flavor oscillations 
and oscillations into sterile states, including the effects of density fluctuations; and 
neutrino-induced and explosive nucleosynthesis. 



1 Introduction 

The lectures summarized here were originally delivered by the authors at a 
summer school held in Canberra, Australia, in January, 1998. The written 
version is a somewhat updated copy of the original, as we have taken this 
opportunity to include some of the exciting neutrino astrophysics results an- 
nounced in summer, 1998. 

The main theme of these lectures is the interplay between the properties of 
the neutrinos — their mass, mixing, and charge conjugation properties — and 
astrophysical phenomena. The first topic is the solar neutrino problem, where 
we describe the current status of the measurements, the physics of the standard 
solar model including helioseismology tests, the input nuclear microphysics, 
and possible solutions. We discuss the MSW mechanism in some detail as well 
as the potential of SuperKamiokande and SNO to demonstrate that oscillations 
are occurring. 

We then discuss Dirac and Majorana neutrino masses, the seesaw mech- 
anism, and possibilities such as pseudoDirac neutrinos. The associated phe- 
nomenology in [3 decay and f3f3 decay is briefly described. 

We discuss a second astrophysical laboratory for neutrino physics, the 
evolution of red giants. Topics include the nuclear physics of the He flash 
and the effects of stellar cooling in red giant and horizontal branch stellar 
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evolution. We illustrated how stellar cooling arguments can place powerful 
constraints on axions and on anomalous properties of the neutrino, such as 
electric and magnetic dipole moments. 

Core-collapse supernovae are the third laboratory. We describe currently 
favored theories of the explosion mechanism and characterize the spectrum 
and flavor of the produced neutrinos. We use supernovae and the solar neu- 
trino problem to motivate a discussion of more exotic aspects of the MSW 
mechanism, including spin-flavor oscillations and the effects of density fluctu- 
ations. We also discuss the nucleosynthesis that can arise because of neutrino 
interactions in the mantle of the star. 

We then turn to the explosive nucleosynthesis associated with the super- 
nova shock wave and the neutrino-driven wind off the proto-neutron star sur- 
face. We discuss the s- and r-processes, and constraints the latter can place 
on neutrino oscillations. These constraints prove relevant to certain terrestrial 
oscillation experiments, such as KARMEN and LSND. 

The audience for these lectures consisted of advanced graduate students 
in nuclear and particle physics: the material is covered at this level and at a 
depth appropriate to a survey. 

2 Solar Neutrinos 

More than three decades ago Ray Davis, Jr. and his collaboratorsi constructed 
a 0.615 kiloton C2CI4 radiochemical solar neutrino detector in the Homestake 
Gold Mine, one mile beneath Lead, South Dakota. Within a few years it was 
apparent that the number of neutrinos detected was considerably below the 
predictions of the standard solar model, that is, the standard theory of main 
sequence stellar evolution. 

Today the results from the CI detector, which have become quite accu- 
rate due to 30 years of careful measurement, have been.-augmented by results 
from four other experiments, the SAGEB and GALLEXcl gallium experiments 
and the Kamiokandetl and SuperKamiokande H water Cerenkov detectors. It 
now appears that the combined results are very difficult to explain — some 
have argued impossible — by any plausible change in the standard solar model 
(SSM). Thus most believe that the answer to the solar neutrino problem is new 
particle physics, most likely some effect like solar neutrino oscillations associ- 
ated with massive neutrinos. With the recent news that SuperKamiokande sees 
direct evidence for oscillations in the azimuthal dependence of atmosphericQ 
neutrinos, it seems that we may be on the threshold of a major discovery. 

The purpose of this first (and longest) lecture is to summarize the solar 
neutrino problem and to present arguments that it represents new particle 
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physics. 



2.1 The Standard Solar Model% 

Solar models trace the evolution of the Sun over the past 4.6 billion years of 
main sequence burning, thereby predicting the present-day temperature and 
composition profiles of the solar core that govern neutrino production. Stan- 
dard solar models share four basic assumptions: 

* The sun evolves in hydrostatic equilibrium, maintaining a local balance 
between the gravitational force and the pressure gradient. To describe this 
condition in detail, one must specify the equation of state as a function of 
temperature, density, and composition. 

* Energy is transported by radiation and convection. While the solar envelope 
is convective, radiative transport dominates in the core region where ther- 
monuclear reactions take place. The opacity depends sensitively on the solar 
composition, particularly the abundances of heavier elements. 

* Thermonuclear reaction chains generate solar energy. The standard model 
predicts that over 98% of this energy is produced from the pp chain conversion 
of four protons into 4 He (see Fig. 1) 

Ap He + 2e+ + 2v c (1) 

with proton burning through the CNO cycle contributing the remaining 2%. 
The Sun is a large but slow reactor: the core temperature, T c ~ 1.5 ■ 10 7 K, 
results in typical center-of-mass energies for reacting particles of ~ 10 keV, 
much less than the Coulomb barriers inhibiting charged particle nuclear re- 
actions. Thus reaction cross sections are small: in most cases, as laboratory 
measurements are only possible at higher energies, cross section data must be 
extrapolated to the solar energies of interest. 

* The model is constrained to produce today's solar radius, mass, and lumi- 
nosity. An important assumption of the standard model is that the Sun was 
highly convective, and therefore uniform in composition, when it first entered 
the main sequence. It is furthermore assumed that the surface abundances of 
metals (nuclei with A > 5) were undisturbed by the subsequent evolution, and 
thus provide a record of the initial solar metallicity. The remaining param- 
eter is the initial He/H ratio, which is adjusted until the model reproduces 
the present solar luminosity after 4.6 billion years of evolution. The resulting 
4 Hc/H mass fraction ratio is typically 0.27 ± 0.01, which can be compared 
to the big-bang value of 0.23 ± 0.01. Note that the Sun was formed from 
previously processed material. 
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Figure 1: The solar pp chain. 

The model that emerges is an evolving Sun. As the core's chemical com- 
position changes, the opacity and core temperature rise, producing a 44% 
luminosity increase since the onset of the main sequence. The temperature 
rise governs the competition between the three cycles of the pp chain: the 
ppl cycle dominates below about 1.6 TO 7 K; the ppll cycle between (1.7-2.3) 
T0 7 K; and the ppIII above 2.4 -10 7 K. The central core temperature of today's 
SSM is about 1.55 -10 7 K. 

The competition between the cycles determines the pattern of neutrino 
fluxes. Thus one consequence of the thermal evolution of our sun is that the 
8 B neutrino flux, the most temperature-dependent component, proves to be 
of relatively recent origin: the predicted flux increases exponentially with a 
doubling period of about 0.9 billion years. 

A final aspect of SSM evolution is the formation of composition gradients 
on nuclear burning timescales. Clearly there is a gradual enrichment of the 
solar core in 4 He, the ashes of the pp chain. Another element, 3 He, is a sort 
of catalyst for the pp chain, being produced and then consumed, and thus 
eventually reaching some equilibrium abundance. The timescale for equilib- 
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rium to be established as well as the eventual equilibrium abundance are both 
sharply decreasing functions of temperature, and thus increasing functions of 
the distance from the center of the core. Thus a steep 3 He density gradient is 
established over time. 

The SSM has had some notable successes. From helioseismology the sound 
speed profile c(r) has been very accurately determined for the outer 90% of the 
Sun, and is in excellent agreement with the SSM. Such studies verify important 
predictions of the SSM, such as the depth of the convective zone. However the 
SSM is not a complete model in that it does not explain all features of solar 
structure, such as the depletion of surface Li by two orders of magnitude. This 
is usually attributed to convective processes that operated at some epoch in 
our sun's history, dredging Li to a depth where burning takes place. 

The principal neutrino-producing reactions of the pp chain and CNO cycle 
are summarized in Table 1. The first six reactions produce j3 decay neutrino 
spectra having allowed shapes with endpoints given by E'" ax . Deviations from 
an allowed spectrum occur for 8 B neutrinos because the 8 Be final state is 
a broad resonance. The last two reactions produce line sources of electron 
capture neutrinos, with widths ~ 2 keV characteristic of the temperature of 
the solar core. Measurements of the pp, 7 Be, and 8 B neutrino fluxes will 
determine the relative contributions of the ppl, ppll, and ppIII cycles to solar 
energy generation. As discussed above, and as later illustrations will show 
more clearly, this competition is governed in large classes of solar models by 
a single parameter, the central temperature T c . The_flux predictions of the 
1998 calculations of Bahcall, Basu, and Pinsonneault (BP98) and of Brun, 
Turck-Chieze and Morel □ arc included in Table 1. 

2.2 Solar Neutrino Detection 

Let us start with a brief reminder about low energy neutrino-nucleus interac- 
tions in detectors. Consider the charged current reaction 

v e + {A,Z)^e- + (A,Z+l) (2) 

Because the momentum transfer to the nucleus is very small for solar neutrinos, 
it can be neglected in the weak propagator, leading to an effective contact 
current-current interaction. If we begin with the simplest case of the free 
neutron decay n — > p + e~ +v e , the corresponding transition amplitude is then 

G F 

Sfi = -^=cos6» c u(p)7 AI (l - g A j 5 )u(n)u(e)-f^{l - ^)u{v) (3) 

where Gf is the weak coupling constant measured in muon decay and cos# c 
gives the amplitude for the weak interaction to connect the u quark to its first- 
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Table 1: Solar neutrino sources and the flux predictions of the Bahcall/Pinsonneault (BP98) 
and Brun/Turck-Chieze/Morel (BTCM98) SSMs in cm"^" 1 . 



Source 




BP98 


BTCM98 


P + P 


^ 2 H + c+ + v 


0.42 


5.94E10 


5.98E10 


13 N- 


13 C + c+ + v 


1.20 


6.05E8 


4.66E8 


15Q _ 


15 N + c+ + v 


1.73 


5.32E8 


3.97E8 




17 + c+ + v 


1.74 


6.33E6 




8 B -» 


s Bo + c+ + v 


- 15 


5.15E6 


4.82E6 


3 He + p 


-^ 4 Hc + c+ + v 


18.77 


2.10E3 




7 Bc + 


c~ -^ 7 Li + v 


0.86 (90%) 


4.80E9 


4.70E9 






0.38 (10%) 






p + e~ 


+ p -^ 2 H + v 


1.44 


1.39E8 


1.41E8 



generation partner, the d quark. The origin of this effective amplitude is the 
underlying standard model predictions for the elementary quark and lepton 
currents. The weak interactions at this level are predicted by the standard 
model to be exactly left handed. Experiment shows that the effective coupling 
of the W boson to the nucleon is governed by 7^(1 — 5.475), as noted above, 
where gA ~ 1.26. The axial coupling is thus shifted from its underlying value 
by the strong interactions responsible for the binding of the quarks within the 
nucleon. 

If an isolated nucleon were the target, one could proceed to calculate the 
cross section from the effective nucleon current given above. The extension 
to nuclear systems traditionally begins with the observation that nucleons in 
the nucleus are rather non-relativistic, v/c <~ 0.1. The amplitude u(p)7 A1 (l — 
<M75)w(^) can be expanded in powers of p/M. The leading vector and axial 
operators are readily found to be 

7o : 1 

7 : p/M ~ v/c 
7075 : 3 ■ p/M ~ v/c 
775 : v 

Thus it is the time-like part of the vector current and the space-like part of 



6 



the axial- vector current that survive in the non-relativistic limit .jj 

If such a non-relativistic reduction is done for our single current one obtains 

S fi ~ cos^(0t(p)0(r^( e ) 7 °(l - lb )u(v) 

-<t> ] '(p)9Ad(t>(ri) ■ u(e)7(l - 75)11(1;)) (4) 

where the 0's are now two-component Pauli spinors for the nucleons. The 
above result can be generalized to include v e reactions by introducing the 
isospin operators t± where r + | n) = | p) and r_ | p) = | n) , with all other 
matrix elements being zero. Thus we can generalize our n — » p amplitude to 
n <-> p by 

^(p)0(n)^^(N)T ± 0(N) 

4>\p)a(j){n) -> cf>i(N)aT±(f)(N). 

This result easily generalizes to nuclear decay. Given our comments about 
exchange currents, the first step is the replacement 

A 

i=l 

crr ± -> y^cr(z)T±(i). 
i=l 

Plugging Sy-i into the standard cross section formula (which involves an average 
over initial and sum over final nuclear spins of the square of the transition 
amplitude) then yields the allowed nuclear matrix element 



2 j^T fl(/llE T ±«H l )l 2 + 5il(/llE ff W T ±«IN)l 2 ) ■ ( 5 ) 



1=1 1=1 



a In a nucleus these currents must be corrected for the presence of meson exchange contri- 
butions. The corrections to the vector charge and axial three-current, which we just pointed 
out survive in the non-relativistic limit, are of order {v/cf ~ 1%. Thus the naive one- 
body currents are a very good approximation to the nuclear currents. In contrast, exchange 
current corrections to the axial charge and vector three-current operators are of order v/c, 
and thus of relative order 1. This difficulty for the vector three-current can be largely cir- 
cumvented, because current conservation as embodied in the generalized Siegert's theorem 
allows one to rewrite important parts of this operator in terms of the vector charge oper- 
ator. In the long-wavelength limit appropriate to solar neutrinos, all terms unconstrained 
by current conservation do not survive. In effect, one has replaced a current operator with 
large two-body corrections by a charge operator with only small corrections. In contrast, the 
axial charge operator is significantly altered by exchange currents even for long-wavelength 
processes like /3 decay. Typical axial-charge /3 decay rates are enhanced by ~ 2 because of 
exchange currents. 
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Our initial calculation for the nucleon treated that particle as structureless. 
Implicitly we assumed that the momentum transfer is much smaller than the 
inverse nucleon size. If we take 10 MeV as a typical solar neutrino momentum 
transfer, these quantities would be in the ratio 1:20. For a light nucleus, the 
corresponding result might be 1:10. This long- wavelength approximation in 
combination with the non-relativistic approximation yields the allowed result, 
where only Fermi and Gamow- Teller operators survive. These are the spin- 
independent and spin-dependent operators appearing above. 

The Fermi operator is proportional to the isospin raising/lowering opera- 
tor: in the limit of good isospin, which typically is good to 5% or better in the 
description of low-lying nuclear states, it can only connect states in the same 
isospin multiplet, that is, states with a common spin-spatial structure. If the 
initial state has isospin (Tj, Mti), this final state has (Tj, Mt% ± 1) for (y, e~) 
and (&, e + ) reactions, respectively, and is called the isospin analog state (IAS). 
In the limit of good isospin the sum rule for this operator in then particularly 
simple 

A A 

£ 2jrfT K/ll ^ T+(i)ll * >|2 = 2jTfT l</ ^ ll ^ T+wl|i>|2 = lN ~ zl (6) 

/ 1 4=1 1 8=1 

The excitation energy of the IAS relative to the parent ground state can be 
estimated accurately from the Coulomb energy difference 

*"~(u5ras- 1J »') Mi,v - (7) 

The angular distribution of the outgoing electron for a pure Fermi (N, Z) + v — > 
(N — 1, Z + 1) + e~ transition is 1 + j3 cos #„ e , and thus forward peaked. Here 
P is the electron velocity. 

The Gamow- Teller (GT) response is more complicated, as the operator 
can connect the ground state to many states in the final nucleus. In general 
we do not have a precise probe of the nuclear GT response apart from weak 
interactions themselves. However a good approximate probe is provided by 
forward-angle (p,n) scattering off nuclei, a technique that has been developed 
in particular by experimentalists at the Indiana University Cyclotron Facilit^O. 
The (p,n) reaction transfers isospin and thus is superficially like (y, e~). At 
forward angles (p,n) reactions involve negligible three-momentum transfers to 
the nucleus. Thus the nucleus should not be radially excited. It thus seems 
quite plausible that forward-angle (p,n) reactions probe the isospin and spin of 
the nucleus, the macroscopic quantum numbers, and thus the Fermi and GT 
responses. For typical transitions, the correspondence between (p,n) and the 
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weak GT operators is believed to be accurate to about 10%. Of course, in a 
specific transition, much larger discrepancies can arise. 

The (p,n) studies demonstrate that the GT strength tends to concentrate 
in a broad resonance centered at a position S = Eqt — Eias relative to the 
IAS given by 

5 ~ ^7.0-28.9 ^ MeV. (8) 

Thus while the peak of the GT resonance is substantially above the IAS for 
N ~ Z nuclei, it drops with increasing neutron excess. Thus 8 ~ for Pb. A 
typical value for the full width at half maximum r is ~ 5 MeV. 

The approximate Ikeda sum rule constrains the difference in the (y, e~) 
and (D,e + ) strengths 

£(|M&(i/, e -)| a " \Mgr(D, e+)| 2 ) = 3(N Z) (9) 

where 

\M{t T (v,e-)\ 2 = ^j^\{f\\J2<rWT + (m\ 2 - (10) 

1 i—l 

In many cases of interest in heavy nuclei, the strength in the (p, e + ) direction 
is largely blocked. For example, in a naive 2sld shell model description of 37 CI, 
the p — ► n direction is blocked by the closed neutron shell at N=20. Thus this 
relation can provide an estimate of the total (3~ strength. Experiment shows 
that the f3~ strength found in and below the GT resonance does not saturate 
the Ikeda sum rule, typically accounting for ~ (60 — 70) % of the total. Mea- 
sured and shell model predictions of individual GT transition strengths tend 
to differ systematically by about the same factor. Presumably the missing 
strength is spread over a broad interval of energies above the GT resonance. 
This is not unexpected if one keeps in mind that the shell model is an approxi- 
mate effective theory designed to describe the long wavelength modes of nuclei: 
such a model should require effective operators, renormalizcdJrom their bare 
values. Phenomenologically, the shell model seems to required ~ 1.0 as 
well as a small spin-tensor term (a £g> Y2(f))j-i of relative strength ~ 0.1. 

The angular distribution of GT (N, Z) + v e -> (N-l, Z+l) + e~ reactions 
is 3 — /3cos#„ e , corresponding to a gentle peaking in the backward direction. 

The above discussion of allowed responses can be repeated for neutral cur- 
rent processes such as (y, The analog of the Fermi operator contributes 
only to elastic processes, where the standard model nuclear weak charge is 
approximately the neutron number. As this operator does not generate tran- 
sitions, it is not yet of much interest for solar or supernova neutrino detection, 
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though there are efforts to develop low-threshold detectors (e.g., cryogenic 
technologies) for recording the modest nuclear recoil energies. The analog of 
the GT response involves 

\M^vy)\ 2 =^j^\(f\\J2<r(^\\i)\ 2 - (11) 

The operator appearing in this expression is familiar from magnetic moments 
and magnetic transitions, where the large isovector magnetic moment ([i v ~ 
4.706) often leads to it dominating the orbital and isoscalar spin operators. 

Finally, there is one purely leptonic reaction of great interest, since it is the 
reaction exploited by Kamiokande and SuperKamiokande. Electron neutrinos 
can scatter off electrons via both charged and neutral current reactions. The 
cross section calculation is straightforward and will not be repeated here. Two 
features of the result are of importance for our later discussions, however. 
Because of the neutral current contribution, heavy-flavor (y^ and v T ) also 
scatter off electrons, but with a cross section reduced by about a factor of 
seven at low energies. Second, for neutrino energies well above the electron 
rest mass, the scattering is sharply forward peaked. Thus this reaction allows 
one to exploit the position of the Sun in separating the solar neutrino signal 
from a large but isotropic background. 

As we mentioned earlier, the first experiment performed was one exploiting 
the reaction 

37 Cl(^,e-) 37 Ar. 

As the threshold for this reaction is 0.814 MeV, the important neutrino sources 
are the 7 Be and 8 B reactions. The 7 Be neutrinos excite just the GT transition 
to the ground state, the strength of which is known from the electron capture 
lifetime of 37 Ar. The 8 B neutrinos can excite all bound states in 37 Ar, including 
the dominant transition to the IAS residing at an excitation of 4.99 MeV. The 
strength of excite-state GT transitions can be determined from the (3 decay 
37 Ca(/3 + ) 37 K, which is the isospin mirror reaction to 37 C1(^, e _ ) 37 Ar. The 
net result is that, for SSM fluxes, 78% of the capture rate should be due 
to 8 B neutrinos, and 15% to 7 Be neutrinos. The measured capture rate l3 
2.56 ±0.16 ± 0.16 SNU (1 SNU = 10~ 36 capture/atom/sec) is about 1/3 the 
standard model value. 

Similar radiochemical experiments were done by the SAGE and GALLEX 
collaborations using a different target, 71 Ga. The special properties of this tar- 
get include its low threshold and an unusually strong transition to the ground 
state of 71 Ge, leading to a large pp neutrino cross section (see Fig. 2). The 
experimental capture rates are 66 ± 13 ± 6 and 76 ± 8 SNU for the SAGE and 
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GALLEX detectors, respectively. The SSM prediction is about 130 SNU0. 
Most important, since the pp flux is directly constrained by the solar lumi- 
nosity in all steady-state models, there is a minimum theoretical value for the 
capture rate of 79 SNU, given standard model weak interaction physics. Note 
there are substantial uncertainties in the 71 Ga cross section due to 7 Be neutrino 
capture to two excited states of unknown strength. These uncertainties were 
greatly reduced by direct calibrations of both detectors using 51 Cr neutrino 
sources. 
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Figure 2: Level scheme for 71 Ge showing the excited states that contribute to absorption of 
pp, 7 Be, 51 Cr, and 8 B neutrinos. 

The remaining experiments, Kamiokande II /III and SuperKamiokande, 
exploited water Cerenkov detectors to view solar neutrinos on an event-by- 
event basis. Solar neutrinos scatter off electrons, with the recoiling electrons 
producing the Cerenkov radiation that is then recorded in surrounding photo- 
tubes. Thresholds are determined by background rates; SuperKamiokande is 



11 



currently operating with a trigger at approximately 6 MeV. The initial ex- 
periment, Kamiokande II/III, found a flux of 8 B neutrinos of (2.91 ±0.24 ± 
0.35) • 10 6 /cm 2 s after about a decade of measurement. Its much larger succes- 
sor SuperKamiokande, with a 22.5 kiloton fiducial volume, yielded the. result 
(2.44 ± 0.05±££ 7 ) • 10 6 /cm 2 s after the first 504 days of measurement^ This 
is about half of the SSM flux. This result continues to improve in accuracy. 



2.3 Uncertainties in Standard Solar Model Parameters 

The pattern of solar neutrino fluxes that has emerged from these experiments 
is 

0(pp)~O.9 SSM ( W ) 
</>( 7 Be) ~ 

0( 8 B) ~ O.43 SSM ( 8 B). (12) 

A reduced 8 B neutrino flux can be produced by lowering the central temper- 
ature of the sun somewhat, as </>( 8 B)~ T c 18 . However, such an adjustment, 
either by varying the parameters of the SSM or by adopting some nonstan- 
dard physics, tends to push the 0( 7 Be)/(/>( 8 B) ratio to higher values rather 
than the low one of eq. (12), 

"( 7 Be) . ^ 10 



? B) 



T- w . (13) 



Thus the observations seem difficult to reconcile with plausible solar model 
variations: one observable, 0( 8 B), requires a cooler core while a second, the 
ratio c/)( 7 Be)/0( 8 B), requires a hotter one. 

An initial question is whether this problem remains significant when one 
takes into account known uncertainties in the parameters of the SSM. While a 
detailed summary of the SSM uncertainties would take us well beyond the lim- 
its of these lectures, a qualitative discussion of pp chain nuclear uncertainties 
is appropriate. This nuclear microphysics has been the focus of a great deal 
of experimental work. The pp chain involves a series of non-resonant charged- 
particle reactions occurring at center-of-mass energies that are well below the 
height of the inhibiting Coulomb barriers. As the resulting small cross sections 
generally preclude laboratory measurements at the relevant energies, one must 
extrapolate higher energy measurements to threshold to obtain solar cross 
sections. This extrapolation is often discussed in terms of the astrophysical 
S-factor 

a(E) = ^-exp(-27T V ) (14) 
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where 77 = ZlZ2a ; with a the fine structure constant and (3 = v/c the relative 
velocity of the colliding particles. This parameterization removes the gross 
Coulomb effects associated with the s-wave interactions of charged, point-like 
particles. The remaining energy dependence of S(E) is gentle and can be ex- 
pressed as a low-order polynomial in E. Usually the variation of S(E) with E is 
taken from a direct reaction model and then used to extrapolate higher energy 
measurements to threshold. The model accounts for finite nuclear size effects, 
strong interaction effects, contributions from other partial waves, etc. As labo- 
ratory measurements are made with atomic nuclei while conditions in the solar 
core guarantee the complete ionization of light nuclei, additional corrections 
must be made to account for the different electronic screening environments. 

Recently a large working group met at a workshop sponsored by the In- 
stitute for Nuclear Theory, University of Washington, to review past work on 
the nuclear reactions of the pp chain and CNO cycle, to recommend best val- 
ues and appropriate errors, and to identify specific issues in experiment and 
theory where additional work is needed. The results are published in Reviews 
of Modern PhysicsL£l. We will not attempt a summary here, but will give one 
or two highlights. 

The most significant recommend change involves-the reaction 7 Be(p, 7) 8 B, 
where the standard Si7(0)~ 22.4 eVb is that givenlill by Johnson et al. Mea- 
surements of Si7(E) are complicated by the need to use radioactive targets and 
thus to determine the areal density of the 7 Be target nuclei. Two techniques 
have been employed, measuring the rate of 478 keV photons from 7 Be decay or 
counting the daughter 7 Li nuclei via the reaction 7 Li (d,p) 8 Li. The low-energy 
data sets for Si7(E) disagree by 25%. This is a systematic normalization prob- 
lem as each data set is consistent with theory in its dependence on E. The 
energy dependence below ~ 500 keV is believed to be quite simple as it is 
determined by the asymptotic nuclear wave function. 

The Seattle working groujxon Si7(E) found that only one low-energy data 
set, that of Filippone et al. Ea, was described in the published literature in 
sufficient detail to be evaluated. The target activity in that experiment had 
been measured by both 478 keV gamma rays and by the (d,p) reaction, with 
consistent results. The resulting recommended value was thus based on this 
measurement, yielding 

5i 7 (0) = 19±^eV b, la. (15) 

A recent measurement is consistent with this value. 

The 3 He(a, 7) 7 Be reaction has been measured by two techniques, by count- 
ing the capture 7 rays and by detecting the resulting 7 Be activity. While the 
two techniques have been used by several groups and have yielded separately 
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consistent results, the capture 7 ray value Si7(0) = 0.507 ± 0.016 keV b is not 
in good agreement with the 7 Be activity value 0.572 ± 0.026 keV-b. The Seat- 
tle working group concluded that the evidence for a systematic discrepancy of 
unknown origin was reasonably strong and recommended that standard proce- 
dures be used in assigning a suitably expanded error. The recommended value 
S 34 (0) is 0.53 ± 0.05. 

These and other recommended values were recently incorporated into the 
BP98 and BTCM98 solar model calculations. While the workshop's recom- 
mended values involve no qualitative changes, there is some broadening of 
error bars. The downward shift in Si7(0) leads to a lower 8 B flux. The work- 
shop's Reviews of Modern Physics article summarizes a substantial amount of 
work on topics not discussed here: screening effects, weak radiative corrections 
to and exchange current effects on p+p, the atomic physics of 7 Be + e~, etc. 
Much of this discussion was useful in evaluating possible uncertainties in solar 
microphysics, and in identifying opportunities for reducing these uncertainties. 

Are uncertainties in the parameters of the SSM a significant source of 
uncertainty? The S-factors discussed above comprise one set of parameters, 
but there are others: the solar lifetime, the opacities, the solar luminosity, etc. 
In order to answer this question while also taking into account correlations 
among the fluxes when input parameters are varied, first Bahcall and UlrichEj 
and later Bahcall and HaxtonEil constructed 1000 SSMs by randomly varying 
five input parameters, the primordial heavy-element-to-hydrogen ratio Z/X 
and S(0) for the p-p, 3 He- 3 He, 3 He- 4 He, and p- 7 Be reactions, assuming for each 
parameter a normal distribution with the mean and standard deviation. (These 
were the parameters assigned the largest uncertainties.) Smaller uncertainties 
from radiative opacities, the solar luminosity, and the solar age were folded 
into the results of the model calculations perturbatively. 

The resulting pattern of 7 Be and 8 B flux predictions is shown in Fig. 3. 
The elongated error ellipses indicate that the fluxes are strongly correlated. 
Those variations producing 0( 8 B) below 0.8</> SSM ( 8 B) tend to produce a re- 
duced 0( 7 Be), but the reduction is always less than 0.8. Thus a greatly reduced 
</>( 7 Be) cannot be achieved within the uncertainties assigned to parameters in 
the SSM. 

A similar exploration, but including parameter variations very far from 
their preferred values, was carried out by Castellani et al. E3, who displayed 
their results as a function of the resulting core temperature T c . The pattern 
that emerges is striking (see Fig. 4) : parameter variations producing the same 
value of T c produce remarkably similar fluxes. Thus T c provides an excellent 
one-parameter description of standard model perturbations. Figure 4 also 
illustrates the difficulty of producing a low ratio of 4>( 7 Be)/4>( 8 B) when T c is 
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Figure 3: SSM 7 Be and S B flux predictions. The dots represent the results of SSM calcu- 
lations where the input parameters were varied according to their assigned uncertainties, as 
described in the text. The 90% and 99% confidence level error ellipses are shown. 



reduced. 

The 1000-solar-model variations were made under the constraint of repro- 
ducing the solar luminosity. Those variations show a similar strong correlation 
with T c 

-Xf^aT- 1 - 2 ^Be)«T c 8 <H 8 B)cxT c 18 . (16) 

Figures 3 and 4 offer a strong argument that reasonable variations in the pa- 
rameters of the SSM, or nonstandard changes in quantities like the metallicity, 
opacities, or solar age, cannot produce the pattern of fluxes deduced from 
experiment (eq. (12)). This would seem to limit possible solutions to errors 
either in the underlying physics of the SSM or in our understanding of neutrino 
properties. 
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Figure 4: The responses of the pp, 'Be, and 8 B neutrino fluxes to the indicated variations in 
solar model input parameters, displayed as a function of the resulting central temperature 

T c . From Castcllani et al. 



2.4 Nonstandard Solar Models 

Nonstandard solar models include both variations of SSM parameters far out- 
side the ranges that are generally believed to be reasonable (some examples 
of which are given in Figure 4), and changes in the underlying physics of the 
model. The solar neutrino problem has been a major stimulus to models: in 
fact, most suggestions were motivated by the hope of producing a cooler Sun 
(T c ~ 0.95T C ) that would avoid conflict with the results of the 37 C1 experi- 
ment. The suggestions included models with low heavy element abundances 
("low Z" models), in which one abandons the SSM assumption that the initial 
heavy element abundances are those we measure today at the Sun's surface; 
periodically mixed solar cores; models where hydrogen is continually mixed 
into the core by turbulent diffusion or by convective mixing; and models where 
the solar core is partially supported by a strong central magnetic field or by its 
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rapid rotation, thereby relaxing the SSM assumption that hydrostatic equilib- 
rium is achieved only through the gas pressure gradient. A larger list is given 
by Bahcall and Davis c3. To illustrate the kinds of consequences such models 
have, two of these suggestions are discussed in more detail below. 

In low-Z models one postulates a reduction in the core metallicity from Z 
~ 0.02 to Z ~ 0.002. This lowers the core opacity (primarily because metals 
are very important to free-bound electron transitions), thus reducing T c and 
weakening the ppll and ppIII cycles. The attractiveness of low-Z models is due 
in part to the existence of additional mechanisms for adding heavier elements 
to the Sun's surface. These include the infall of comets and other debris, as 
well as the accumulation of dust as the Sun passes through interstellar clouds. 
However, the increased radiative energy transport in low-Z models leads to a 
thin convective envelope, in contradiction to interpretations of the 5-minute 
solar surface oscillations. A low He mass fraction also results. As diffusion of 
material from a thin convective envelope into the interior would deplete heavy 
elements at the surface, investigators have also questioned whether present 
abundances could have accumulated in low-Z models. Finally, the general 
consistency of solar heavy element abundances with those observed in other 
main sequence stars makes the model appear contrived. 

Models in which the solar core (~ 0.2 M Q ) is intermittently mixed break 
the standard model assumption of a steady-state Sun: for a period of sev- 
eral million years (the thermal relaxation time for the core) following mixing, 
the usual relationship between the observed surface luminosity and rate of 
energy (and neutrino) production is altered as the Sun burns out of equilib- 
rium. Calculations show that both the luminosity and the 8 B neutrino flux 
are suppressed while the Sun relaxes back to the steady state. Such models 
have been considered seriously because of instabilities associated with large 
gradients in the 3 He abundance, which in equilibrium varies as ~ T~ 6 , where 
T is the local temperature. The resulting steep profile is unstable under finite 
amplitude displacements of a volume to smaller r: the energy released by the 
increased 3 He burning at higher T can exceed the energy in the perturbation. 
For a discussion of the plausibility of such_a trigger for core mixing, one can 
see the original work of Dilke and Gough Ell as well as a more recent critique 
by Merryfieldu. The possibility that continuous mixing on time scales of 3 Hc 
mixing could produce a flux pattern close to that observed (e.g., a suppression 
in both the 8 B neutrino flux and the 7 Be/ 8 B flux ratio) was recently discussed 
by Cumming and Haxtone3. 

This discussion of two of the more seriously explored nonstandard model 
possibilities illustrates how changes motivated by the solar neutrino problem 
often produce other, unwanted consequences. In particular, many experts feel 



17 



that the good SSM agreement with helioseismology is likely to be destroyed 
by changes such as those discussed above. 

Figure 5 is an illustration by Hata et al.EH of the flux predictions of several 
nonstandard models, including a low-Z model consistent with the 37 C1 results. 
As in the Castellani et al. exploration, the results cluster along a track that 
defines the naive T c dependence of the 0( 7 Be)/</>( 8 B) ratio, well separated from 
the experimental contours. 




(K 7 Be)/(M 7 Be) SSM 

Figure 5: The fluxes allowed by the combined results of the various solar neutrino experi- 
ments compared to the results of SSM variations and various nonstandard solar models. The 
solid line in the naive T c power law discussed in the text. From Hata et al. 

There is now a popular argument that no such nonstandard model can 
solve the solar neutrino problem: if one assumes undistorted neutrino spectra, 
no combination of pp, 7 Be, and 8 B neutrino fluxes fits the experimental results 
wello. In fact, in an unconstrained fit, the required 7 Be flux is unphysical, 
negative by almost 3<t. Thus, barring some unfortunate experimental error, it 
appears we are forced to look elsewhere for a solution. 
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If experimental error, SSM parameter uncertainties, and nonstandard solar 
physics are ruled out as potential solutions, new particle physics is left as the 
leading possibility. Suggested particle physics solutions of the solar neutrino 
problem include neutrino oscillations, neutrino decay, neutrino magnetic mo- 
ments, and weakly interacting massive particles. Among these, the Mikhcycv- 
Smirnov-Wolfenstein effect — neutrino oscillations enhanced by matter inter- 
actions — is widely regarded as the most plausible. 



2. 5 Helio seismology 

Earlier it was mentioned that measurements of the sound velocity within the 
Sun, deduced from observations of surface oscillations, provide a powerful check 
on the SSM. In this section the basic physics of helioseismology is reviewed. 

A static, stable star at spherically-symmetric equilibrium can be charac- 
terized with pressure p(r), mass density p(r), the gravitational potential 4>(r), 
the rate of nuclear energy generation e(r), temperature T(r), the energy flux 
F and the entropy s. Introducing the adiabatic indices 



and the total derivative 



Dt dt 

one can write down the equation of motion 



<91ogp 

9 log 
d log p 



° a +W (19) 



r.-i-r3?T|. <»> 



Dw 

p— = - Vp - pV0, (20) 



the equation of continuity 



^ + P V-v = 0, (21) 
Poisson's equation for gravitational attraction 

V 2 = AitGp, (22) 
and an equation describing energy conservation 

I^- ri i^ = ^(pe-V.P). (23) 
p Dt p Dt p VP ; v ; 

19 



These equations describe a static star. To do stellar seismology one introduces 
Eulerian (i.e. at a given point) perturbations on the physical quantities, e.g. 

p{r,t)=p (r)+p'(r,t) (24) 

where po denotes the equilibrium value and the displacement is calculated from 
the velocity amplitude 

v = |(^r). (25) 

Inserting expressions like Eq. (^4|) in the above equations and subtracting the 
equilibrium equations one obtains equations that describe the perturbations. 
From the conservation of momentum one gets 



&Sr = _ v j + £ v v<// ^ (26) 
at* p 



The equation of continuity gives 

p' + V • (p6r) = 0. (27) 

Poisson's equation becomes 

VV = AnGp' , (28) 

and the energy equation yields 

^ + -Sr ■ Vp = -I- (- + -Sr ■ VpV (29) 

pp ri \p p j 

In these equations for convenience we dropped the subscript zero in writing 
down the equilibrium values. To obtain the normal modes of a star we assume 
a time dependence of exp(— iut) for the perturbations: 

p'(r,t) ~ p'(r)Y em (6,<f>)ex.p(-iLut). (30) 

Using Eq. (|30| ) and introducing the auxiliary quantity 

#(r) = c 2 p 1/2 V • Sr, (31) 

where c is the adiabatic sound speed 

c = , (32) 

9 
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Eqs. (26) through (29) can be written in a compact form 



d 2 ^ 
dr 2 



+ 




CO 



,2 



£{£+ l)c 2 




* ~ 0. 



(33) 



In Eq. ([33]) we used the buoyancy frequency: 

^2 = Gm (r) f 1 dlogp _ dlogp 
r \Yi dr dr 



(34) 



and the acoustical cut-off frequency 



CO 



,2 




(35) 



where the density scale height is 



H = -(dlogp/dr)- 1 



(36) 



One observes that oscillation frequencies are determined by the sound speed 
profile. The sound speed profile calculated using the Bahcall-Pinsonneault 
1998 SSM is shown in Figure 6. 

Eq. ( |33l ) resembles the Schroedinger equation in quantum mechanics. For 
the present Sun the buoyancy frequency is approximately constant in the ra- 
diative zone except very near the core, but vanishes in the convective zone. 
The acoustical cut-off frequency is a monotonically decreasing function of 
the distance from the center of the Sun. Inserting this "potential" to the 
"Schroedinger equation", Eq. (|33|), one observes that for N 2 /ui 2 <C 1 the am- 
plitude of the oscillations die out in the radiative zone ( "classically forbidden" ) 
of the Sun. The resulting oscillations are confined to the convective (outer) 

zone of the Sun and called p- modes (for pressure). For e ^t 1 \ c ^> 1 the situ- 
ation is reversed. The "classically allowed" region is the core of the Sun and 
the amplitudes die out in the convective zone. These oscillations are called 
g- modes (for gravity). Note that g-mode oscillation amplitudes vanish at the 
solar surface, hence it is very difficult to directly observe g-modes. On the 
other hand p-mode oscillations are readily measurable by observing the solar 
surface. Eq. ( |33| ) indicates that p-mode oscillations with different i values 
penetrate to different depths, the observed frequencies ujg are determined by 
conditions at different parts of the Sun. 

It is possible to gain insight to the properties of solar oscillations by regard- 
ing the equations outlined above as an eigenvalue problem in a linear Hilbert 
space. Hence it is possible to directly relate perturbations in the sound speed to 
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Figure 6: Sound speed profile in the Sun. From Bahcall and Pinsonneault. 



the perturbations in u>£. One starts with the sound speed profile and oscillation 
frequencies calculated in a reference solar model. These quantities are taken as 
the unperturbed quantities. Then using standard Raleigh-Schroedinger per- 
turbation theory one relates the difference between observed and calculated 
frequencies to the deviation of the sound speed from the model prediction. 
A very readable introduction the,the theoretical aspects of helioseismology is 
available on the world wide webcH There is very good agreement between cal- 
culated and observed sound speeds in the Sun. Figure 7 shows the fractional 
difference between the predicted and observed sound speed profiles u. Sound 
speed profiles deduced from helioseismology provides an important constraint 
on solar models. 
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Figure 7: Fractional difference between the predicted and observed sound speed profiles. 
From Bahcall, Basu, and Pinsonneault. 



2.6 Neutrino Oscillations 

One odd feature of particle physics is that neutrinos, which are not required by 
any symmetry to be massless, nevertheless must be much lighter than any of 
the other known fcrmions. For instance, the current limit on the v e mass is < 5 
cV. The standard model requires neutrinos to be massless, but the reasons are 
not fundamental. Dirac mass terms raj), analogous to the mass terms for other 
fermions, cannot be constructed because the model contains no right-handed 
neutrino fields. Neutrinos can also have Majorana mass terms 

v c L mLVL and v^tcirvr, (37) 

where the subscripts L and R denote left- and right-handed projections of 
the neutrino field v, and the superscript c denotes charge conjugation. The 
first term above is constructed from left-handed fields, but can only arise as 
a nonrenormalizable effective interaction when one is constrained to generate 
niL with the doublet scalar field of the standard model. The second term is 
absent from the standard model because there are no right-handed neutrino 
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fields. 

None of these standard model arguments carries over to the more general, 
unified theories that theorists believe will supplant the standard model. In 
the enlarged multiplets of extended models it is natural to characterize the 
fermions of a single family, e.g., v e , e, u, d, by the same mass scale mp. Small 
neutrino masses are then frequently explained as a result of the Majorana 
neutrino masses. In the seesaw mechanisnJa, 

M^(° T mD ). (38) 

Diagonalization of the mass matrix produces one light neutrino, mi; g ht ~ 
and one unobservably heavy, mhcavy ~ tur. The factor (mo/ 'ttir) is the needed 
small parameter that accounts for the distinct scale of neutrino masses. The 
masses for the v e ^v^ and v T are then related to the squares of the correspond- 
ing quark masses m u , m c , and m t . Taking m R ~ 10 16 GeV, a typical grand 
unification scale for models built on groups like SO (10), the seesaw mechanism 
gives the crude relation 

m Ve : m Vli : m„ T <-> 2 • 10~ 12 : 2 • 1(T 7 : 3 • l(T 3 eV. (39) 

The fact that solar neutrino experiments can probe small neutrino masses, and 
thus provide insight into possible new mass scales m R that are far beyond the 
reach of direct accelerator measurements, has been an important theme of the 
field. 

One of the most interesting possibilities for solving the solar neutrino prob- 
lem has to do with neutrino masses. For simplicity we will discuss just two 
neutrinos. If a neutrino has a mass m, we mean that as it propagates through 
free space, its energy and momentum are related in the usual way for this mass. 
Thus if we have two neutrinos, we can label those neutrinos according to the 
eigenstates of the free Hamiltonian, that is, as mass eigenstates. 

But neutrinos are produced by the weak interaction. In this case, we have 
another set of eigenstates, the flavor eigenstates. We can define a v e as the 
neutrino that accompanies the positron in j3 decay. Likewise we label by 
the neutrino produced in muon decay. 

Now the question: are the eigenstates of the free Hamiltonian and of the 
weak interaction Hamiltonian identical? Most likely the answer is no: we know 
this is the case with the quarks, since the different families (the analog of the 
mass eigenstates) do interact through the weak interaction. That is, the up 
quark decays not only to the down quark, but also occasionally to the strange 
quark. (This is why we had a cos6* c in our weak interaction amplitude: the 
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amplitude for u — > s is proportional to sin# c .) Thus we suspect that the weak 
interaction and mass eigenstates, while spanning the same two-neutrino space, 
are not coincident: the mass eigenstates \i>i) and 1 1^2) (with masses mi and 
m-i) are related to the weak interaction eigenstates by 

\u e ) = cos6 v \vi) + sin6»„|^ 2 ) 

= -sm0„k) +cos(9„|ia>) (40) 

where V is the (vacuum) mixing angle. 

An immediate consequence is that a state produced as a \u e ) or a |^,) at 
some time t — for example, a neutrino produced in (3 decay — does not remain 
a pure flavor eigenstate as it propagates away from the source. This is because 
the different mass eigenstates comprising the neutrino will accumulate differ- 
ent phases as they propagate downstream, a phenomenon known as vacuum 
oscillations (vacuum because the experiment is done in free space). To see 
the effect, suppose we produce a neutrino in some (3 decay where we measure 
the momentum of the initial nucleus, final nucleus, and positron. Thus the 
outgoing neutrino is a momentum eigenstateEj. At time t—Q 

\v(t = 0)) = \u e ) =cos0,|i<i)+flra0 t ,|i*). (41) 

Each eigenstate subsequently propagates with a phase 

^(k-x—ujt) _ ^i{k-x — ^/m?+k 2 t) (42) 

But if the neutrino mass is small compared to the neutrino momentum/energy, 
one can write 

^ m 2 + k^k(l + ^). (43) 

Thus we conclude 

\v(ij) — e i ^- s - kt -( m2 i+ m l) t / 4k ) 

x[cos0 v \ Vl )e lSm2t/4k + smO v \v 2 )e- lSm2t / 4k }. (44) 

We see there is a common average phase (which has no physical consequence) 
as well as a beat phase that depends on 

5m 2 = m\ — mf. (45) 

Now it is a simple matter to calculate the probability that our neutrino state 
remains a \v e ) at time t 

P Ve {t) = \{v e \v{t))? 

{5m 2 A f 
= f - sin 2 20 v sin 2 — — ^ 1 - - sin 2 29 v (46) 
\ ik J 2 
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where the limit on the right is appropriate for large t. Now E ~ k, where E 
is the neutrino energy, by our assumption that the neutrino masses are small 
compared to k. Thus we can reinsert the units above to write the probability 
in terms of the distance x of the neutrino from its source, 



(When one properly describes the neutrino state as a wave packet, the large- 
distance behavior follows from the eventual separation of the mass eigenstates.) 
If the the oscillation length 



is comparable to or shorter than one astronomical unit, a reduction in the solar 
u e flux would be expected in terrestrial neutrino oscillations. 

The suggestion that the solar neutrino problem could be explained by 
neutrino oscillations was first made by Pontecorvo in 1958, who pointed out 
the analogy with Kq <-> Kq oscillations. From the point of view of particle 
physics, the sun is a marvelous neutrino source. The neutrinos travel a long 
distance and have low energies (~ 1 MeV), implying a sensitivity down to 



In the seesaw mechanism, Sm 2 ~ m 2 ., so neutrino masses as low as to 2 <~ 
1CP 6 eV could be probed. In contrast, terrestrial oscillation experiments with 
accelerator or reactor neutrinos are typically limited to Sm 2 > 0.1 cV 2 . 

From the expressions above one expects vacuum oscillations to affect all 
neutrino species equally, if the oscillation length is small compared to an as- 
tronomical unit. This is somewhat in conflict with the solar neutrino data, as 
we have argued that the 7 Be neutrino flux is quite suppressed. Furthermore, 
there is a weak theoretical prejudice that 9 V should be small, like the Cabibbo 
angle. The first objection, however, can be circumvented in the case of "just 
so" oscillations where the oscillation length is comparable to one astronomical 
unit. In this case the oscillation probability becomes sharply energy depen- 
dent, and one can choose dm 2 to preferentially suppress one component (e.g., 
the monochromatic 7 Be neutrinos). This scenario has been explored by sev- 
eral groups and remains an interesting possibility. However, the requirement 
of large mixing angles remains. 

Below we will see that stars allow us to "get around" this problem with 
small mixing angles. In preparation for this, we first present the results above 




(47) 



dm 2 > 10- 12 eV 2 . 



(49) 
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in a slightly more general way. The analog of Eq. ( [44| ) for an initial muon 
neutrino (\v(t = 0)) = 1^)) is 

\v(t)) = e i (k-£-kt-(ml+ml)t/4k) 

x[-sin#> 1 )e i5m2t / 4 ' £ +cose v \u 2 )e~ lSm2t/ik ] (50) 

Now if we compare Eqs. (|44|) and we see that they are special cases of a 
more general problem. Suppose we write our initial neutrino wave function as 

\u{t = 0)) = a e (t = 0)\v e ) + a^t = 0)|^). (51) 

Then Eqs. (fS]) and (|5(]) tell us that the subsequent propagation is described 
by changes in a e (x) and a^{x) according to (this takes a bit of algebra) 

. d ( a e \ _ 1 f -8m 2 cos29 v 8m 2 sin 29 v \ ( a e \ ,r \ 

1 dx ) ~ 4E \ 8m 2 sin29 v Sm 2 cos 2B V ) \ afl J ' ^ ' 

Note that the common phase has been ignored: it can be absorbed into the 
overall phase of the coefficients a e and a^, and thus has no consequence. Also, 
we have equated x — t, that is, set c = 1. 



2. 7 The Mikheyev-Smirnov- Wolfenstein Mechanism 

The view of neutrino oscillations changed when Mikheyev and Smirnov El 
showed in 1985 that the density dependence of the neutrino effective mass, 
a phenomenon first discussed by Wolfenstein in 1978, could greatly enhance 
oscillation probabilities: a v e is adiabatically transformed into a as it tra- 
verses a critical density within the sun. It became clear that the sun was not 
only an excellent neutrino source, but also a natural regenerator for cleverly 
enhancing the effects of flavor mixing. 

While the original work of Mikheyev and Smirnov was numerical, their 
phenomenon was soon understood analytically as a level-crossing problem. If 
one writes the neutrino wave function in matter as in Eq. (|5l|), the evolution 
of a e (x) and a M (a;) is governed by 

.d_fa e \ _ J_ / 2EV2G F p(x) - 8m 2 cos 29 v Sm 2 sin 29 v \ / a e \ 
l dx \ a^J ~ 4~E \8m 2 sin29 v - 2E^2G F p[x) + 8m 2 cos 29 v ) ) 

(53) 

where Gf is the weak coupling constant and p{x) the solar electron density. 
If p(x) = 0, this is exactly our previous result and can be trivially integrated 
to give the vacuum oscillation solutions of Sec. 2.5. The new contribution 
to the diagonal elements, 2E\[2Gfp{x), represents the effective contribution 
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to Ml that arises from neutrino-electron scattering. The indices of refraction 
of electron and muon neutrinos differ because the former scatter by charged 
and neutral currents, while the latter have only neutral current interactions. 
The difference in the forward scattering amplitudes determines the density- 
dependent splitting of the diagonal elements of the new matter equation. 

ft is helpful to rewrite this equation in a basis consisting of the light and 
heavy local mass eigenstates (i.e., the states that diagonalize the right-hand 
side of the equation) , 

\vl{x)) = cos9(x)\v e ) - sa\9{x)\v^) 

\vh{x)) = sm9(x)\v e ) +cos6(x)\v li ). (54) 
The local mixing angle is defined by 

sin 29 y 



sm29(x) = 



^Jx 2 (x)+sin 2 29 v 



cos2fl(x) = X ^ = (55) 



yjx 2 {x) +sin 2 29 v 



where X{x) = 2\[2G 'f } p{x)E / '8m 2 - cos2# w . Thus 9{x) ranges from 9 V to tt/2 
as the density p(x) goes from to oo. 
If we define 

\v(x)) = a H (x)\v H (x)) + a L (x)\v L (x)), (56) 

the neutrino propagation can be rewritten in terms of the local mass eigenstates 

{ d_ fa H \ = f \{x) ia(x) Wa fl \ ^ 
dx\a L ) \-ia(x) -X(x) J \a L J 

with the splitting of the local mass eigenstates determined by 

2\(x) - S -^Jx 2 (x)+sm 2 29 v (58) 
and with mixing of these eigenstates governed by the density gradient 

= (E_\ V2G F £ P (x) S in29 v 
a[X) \5m 2 ) X 2 (x)+^29 v ■ 159 ' 

The results above are quite interesting: the local mass eigenstates diagonalize 
the matrix if the density is constant. In such a limit, the problem is no more 
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complicated than our original vacuum oscillation case, although our mixing 
angle is changed because of the matter effects. But if the density is not con- 
stant, the mass eigenstates in fact evolve as the density changes. This is the 
crux of the MSW effect. Note that the splitting achieves its minimum value, 
sin 20„, at a critical density p c — p(x c ) 

2V2EG F Pc = Sm 2 cos 20„ (60) 

that defines the point where the diagonal elements of the original flavor matrix 
cross. 

Our local-mass-eigenstate form of the propagation equation can be trivially 
integrated if the splitting of the diagonal elements is large compared to the off- 
diagonal elements (see Eq. (57)), 



X(x) 



a(x) 



sin 2 26^ 5m 2 1 [X(x) 2 + sin 2 20„] 3 / 2 
cos20 w ~2E J±MeT\ sin 3 20 v 

I p c ax I 



» 1, (61) 



a condition that becomes particularly stringent near the crossing point, 

sin 2 20„ Sm 2 1 



7c = l(x c ) 



cos 20„ 2E 



dp(x) I _ 
dx < x — Xc 



> 1. (62) 



The resulting adiabatic electron neutrino survival probability valid when 
7c > 1, is 

P; c diab = i + i cos 20„ cos 20, (63) 

where 9i = 9{xi) is the local mixing angle at the density where the neutrino 
was produced. 

The physical picture behind this derivation is illustrated in Figure 8. One 
makes the usual assumption that, in vacuum, the v e is almost identical to the 
light mass eigenstate, z^l(0), i.e., m\ < 777,2 and cos0„ ~ 1. But as the density 
increases, the matter effects make the v e heavier than the f M , with v e — ► vh{x) 
as p(x) becomes large. The special property of the Sun is that it produces u e s 
at high density that then propagate to the vacuum where they are measured. 
The adiabatic approximation tells us that if initially v & ~ vjj (x) , the neutrino 
will remain on the heavy mass trajectory provided the density changes slowly. 
That is, if the solar density gradient is sufficiently gentle, the neutrino will 
emerge from the sun as the heavy vacuum eigenstate, ~ v^. This guarantees 
nearly complete conversion of is e s into v^s, producing a flux that cannot be 
detected by the Homestake or SAGE/GALLEX detectors. 
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Figure 8: Schematic illustration of the MSW crossing. The dashed lines correspond to the 
electron-electron and muon-muon diagonal elements of the M 2 matrix in the flavor basis. 
Their intersection defines the level-crossing density p c . The solid lines are the trajectories 
of the light and heavy local mass eigenstates. If the electron neutrino is produced at high 
density and propagates adiabatically, it will follow the heavy-mass trajectory, emerging from 

the sun as a 



But this does not explain the curious pattern of partial flux suppressions 
coming from the various solar neutrino experiments. The key to this is the 
behavior when 7 C < 1. Our expression for 7(2;) shows that the critical region 
for non-adiabatic behavior occurs in a narrow region (for small 6 V ) surrounding 
the crossing point, and that this behavior is controlled by the derivative of 
the density. This suggests an analytic strategy for handling non-adiabatic 
crossings: one can replace the true solar density by a simpler (intcgrable!) 
two-parameter form that is constrained to reproduce the true density and its 
derivative at the crossing point x c . Two convenient choices are the linear 
(p(x) = a + bx) and exponential (p(x) = ae~ bx ) profiles. As the density 
derivative at x c governs the non-adiabatic behavior, this procedure should 
provide an accurate description of the hopping probability between the local 
mass eigenstates when the neutrino traverses the crossing point. The initial 
and ending points Xi and x/ for the artificial profile are then chosen so that 
p(xi) is the density where the neutrino was produced in the solar core and 
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p(xf) — (the solar surface), as illustrated in in Figure 9. Since the adiabatic 
result (padiab^ depends only on the local mixing angles at these points, this 
choice builds in that limit. But our original flavor-basis equation can then be 
integrated exactly for linear and exponential profiles, with the results given in 
terms of parabolic cylinder and Whittaker functions, respectively. 
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Figure 9: The top figure illustrates, for one choice of sin 2 29 and 8m 2 , that the region of 
non-adiabatic propagation (solid line) is usually confined to a narrow region around the 
crossing point r c . In the lower figure, the solid lines represent the solar density and a 
linear approximation to that density that has the correct initial and final values, as well 
as the correct density and density derivative at r c . Thus the linear profile is a very good 
approximation to the sun in the vicinity of the crossing point. The MSW equations can 
be solved analytically for this wedge. By extending the wedge to ±oo (dotted lines) and 
assuming adiabatic propagation in these regions of unphysical density, one obtains the simple 
Landau-Zener result discussed in the text. 



That result can be simplified further by observing that the non-adiabatic 
region is generally confined to a narrow region around x c , away from the end- 
points Xi and Xf. We can then extend the artificial profile to x — ±oo, as 
illustrated by the dashed lines in Figure 9. As the neutrino propagates adi- 
abatically in the unphysical region x < Xi, the exact solution in the physical 
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region can be recovered by choosing the initial boundary conditions 



I \ i \ n _i I K x ) dx 

a^(— ooj = — CLu{— oo) = cost^e J -°° 

a#(— oo) = a e (— oo) = sm9ie J -» . (64) 

That is, oo)) will then adiabatically evolve to |z/(xj)) = \v e ) as x goes from 
— oo to Xi. The unphysical region x > x/ can be handled similarly. 

With some algebra a simple generalization of the adiabatic result emerges 
that is valid for all 5m 2 /E and 8 V 

P»e = \ + \ cos 20 v cos 20»(1 - 2P hop ) (65) 

where Phop is the Landau-Zener probability of hopping from the heavy mass 
trajectory to the light trajectory oxl traversing the crossing point. For the 
linear approximation to the densityOtj, 

Pt P = (66) 

As it must by our construction, P Uc reduces to Pj^ lab for j c 3> 1. When the 
crossing becomes non-adiabatic (e.g., 7 C <C 1 ), the hopping probability goes 
to 1, allowing the neutrino to exit the sun on the light mass trajectory as a v ei 
i.e., no conversion occurs. 

Thus there are two conditions for strong conversion of solar neutrinos: 
there must be a level crossing (that is, the solar core density must be sufficient 
to render v e ~ vnixi) when it is first produced) and the crossing must be 
adiabatic. The first condition requires that Sm 2 /E not be too large, and the 
second 7 C > 1. The combination of these two constraints, illustrated in Fig. 

10, defines a triangle of interesting parameters in the ^ sin 2 29 v plane, 

as Mikheyev and Smirnov found by numerically integration. A remarkable 
feature of this triangle is that strong v e — > conversion can occur for very 
small mixing angles (sin 2 26 ~ 10~ 3 ), unlike the vacuum case. 

One can envision superimposing on Fig. 10 the spectrum of solar neutri- 
nos, plotted as a function of ^g- for some choice of 5m 2 . Since Davis sees 
some solar neutrinos, the solutions must correspond to the boundaries of the 

c 2 

triangle in Fig. 10. The horizontal boundary indicates the maximum ^g- 
for which the sun's central density is sufficient to cause a level crossing. If 
a spectrum properly straddles this boundary, we obtain a result consistent 
with the Homestake experiment in which low energy neutrinos (large 1/E) lie 
above the level-crossing boundary (and thus remain ^ e 's), but the high-energy 
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Figure 10: MSW conversion for a neutrino produced at the sun's center. The upper shaded 
region indicates those 8m? /E where the vacuum mass splitting is too great to be overcome 
by the solar density. Thus no level crossing occurs. The lower shaded region defines the 
region where the level crossing is non-adiabatic (7 C less than unity). The unshaded region 
corresponds to adiabatic level crossings where strong u e — ► 



neutrinos (small 1/E) fall within the unshaded region where strong conversion 
takes place. Thus such a solution would mimic nonstandard solar models in 
that only the 8 B neutrino flux would be strongly suppressed. The diagonal 
boundary separates the adiabatic and non-adiabatic regions. If the spectrum 
straddles this boundary, we obtain a second solution in which low energy neu- 
trinos lie within the conversion region, but the high-energy neutrinos (small 
1/E) lie below the conversion region and are characterized by 7 <C 1 at the 
crossing density. (Of course, the boundary is not a sharp one, but is charac- 
terized by the Landau-Zener exponential). Such a non-adiabatic solution is 
quite distinctive since the flux of pp neutrinos, which is strongly constrained 
in the standard solar model and in any steady-state nonstandard model by 
the solar luminosity, would now be sharply reduced. Finally, one can imag- 
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inc "hybrid" solutions where the spectrum straddles both the level-crossing 
(horizontal) boundary and the adiabaticity (diagonal) boundary for small 8, 
thereby reducing the 7 Be neutrino flux more than either the pp or 8 B fluxes. 

What are the results of a careful search for MSW solutions satisfying 
the Homestakc, Kamiokandc/SupcrKamiokande, and SAGE/GALLEX con- 
straints? This has been explored in detail by several groups. One solution, 
corresponding to a region surrounding Sm 2 <~ 6-lCT 6 cV 2 and sin 2 28 v <~ 6-10~ 3 , 
is the hybrid case described above. It is commonly called the small-angle solu- 
tion. A second, large-angle solution exists, corresponding to Sm 2 ~ 10 _5 eV 2 
and sin 2 28 v ^0.6. These solutions can be distinguished by their characteristic 
distortions of the solar neutrino spectrum. The survival probabilities P^ ISW (E) 
for the small- and large-angle parameters given above are shown as a function 
of E in Fig. 11. 
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Figure 11: MSW survival probabilities P(E„) for typical small angle and large angle solutions. 



The MSW mechanism provides a natural explanation for the pattern of 
observed solar neutrino fluxes. While it requires profound new physics, both 
massive neutrinos and neutrino mixing are expected in extended models. The 
small-angle solution corresponds to Sm 2 ~ 10~ 5 eV 2 , and thus is consistent 
with 77i2 ~ few -10~ 3 cV. This is a typical v T mass in models where thr <~ 
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77ZGUT- This mass is also reasonably close to atmospheric neutrino values. 
On the other hand, if it is the participating in the oscillation, this gives 
■niR ~ 10 12 GeV and predicts a heavy v T ~ 10 eV. Such a mass is of great 
interest cosmologically as it would have consequences for supernova physics, 
the dark matter problem, and the formation of large-scale structure. 



2.8 SuperKamiokande, SNO, and the MS W Mechanism 

SuperKamiokande and Sudbury Neutrino Observatory (SNO) detectors are 
real-time counting detectors in contrast to the radiochemical detectors such 
as Homestake, GALLEX, and SAGE, which can only determine a time- and 
energy-integral of the flux (cf. Section 2.2). Both the SuperKamiokande and 
SNO can detect neutrinos through the reaction 

v x + e~ — > v x + e~ . (67) 

The electrons coming from this reaction are confined to a forward cone. Hence 
detecting the Cerenkov radiation from the final electron one can determine 
neutrino's direction. In this reaction it is very difficult to determine the energy 
of the neutrino from the measured energy of the final electron because of the 
kincmatical broadening. However the measured energy spectra of the recoil 
electrons can nevertheless yield valuable information about the neutrino energy 
spectrum. Indeed recent SuperKamiokande measurement of the recoil electron 
energy spectrum is consistently below the SSM predictions at electron energies 
up to 14 MeV (which itself is consistent with other neutrino experiments). 
However for the higher energy bins the data are anomalously high as compared 
to the predictions. Although it can be fitted by assuming a very large hep 
neutrino fluxL3 the origin of this anomalous behavior is not understood. 

In addition to the reaction Eq. ( |67| ) SNO can detect neutrinos through 
the charged-current reaction 

v e + d — > p + p + e~ , (68) 

and the neutral current reaction 

v x (p x ) + d — > v x (px) +P + n. (69) 

The neutrons produced in Eq. ( |69|) can be detected either using (n, 7) reactions 
on salt dissolved in the heavy water or by using 3 He proportional counters. The 



electrons coming from the reaction (68) are essentially monochromatic with 



energies ~ E v — 1.44 MeV and they have a very different angular distribution 
(1 — cos6* e /3) with respect to the neutrino direction than that of the electrons 
coming from the reaction ( |67j ) (which are constrained to the forward cone). 
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Figure 12: The dashed line shows the spectrum distortion at SNO for the small-angle MSW 
solution (8m? ~ 5 X 10 — 6 eV 2 and sin 29 ~ 0.01). The solid line is the spectrum without 
MSW oscillations, normalized to the same total rate as with MSW oscillations. 



It should be emphasized that since matter-enhanced neutrino conversion is 
energy-dependent, MSW mechanism will distort the neutrino energy spectra. 
However since only the energy of the final state electron can be measured, if 
a broad range of electron energies correspond to a given neutrino energy this 
distortion will be smeared. The charged current break-up reaction, Eq. (pq), 
is more suitable for this purpose since the final electron energy is distributed 
around a very narrow peak centered at the initial neutrino energy. Spectrum 
distortion at SNO for the small-angle MSW solution (5m 2 ~ 5 x 10~ 6 cV 2 and 
sin 26 ~ 0.01) is shown in Figure 12. 

The reaction ([39|) measures the total flux of all the neutrinos that par- 
ticipate in the low-energy weak interactions. Hence if there are no sterile 
neutrinos (more about this later) it can be used to test the neutrino oscillation 
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hypothesis. For example if the electron neutrinos are being converted to muon 
neutrinos in the Sun, flux measured using the reaction (|6^) should agree with 
the SSM predictions. Finally it is worth mentioning that the electron antineu- 
trinos coming from a supernova or a possible spin-flavor precession in the Sun 
(cf. Section 2.9) can also be detected through the reaction 

V e + d — >• n + n + e+. (70) 



2.9 Other Particle Physics Solutions 

If the MSW mechanism proves not to be the solution of the solar neutrino prob- 
lem, it still will have greatly enhanced the importance of solar neutrino physics: 
the existing experiments have ruled out large regions in the 5m 2 — sin 2 28 v plane 
(corresponding to nearly complete v e — > Vn conversion) that remain hopelessly 
beyond the reach of accelerator neutrino oscillation experiments. 

A number of other particle physics solutions have been considered, such as 
neutrino decay, matter-catalyzed neutrino decay, and solar energy transport 
by weakly interacting massive particles. But perhaps the most interesting 
possibility, apart from the MSW mechanism, was stimulated by suggestions 
that the 37 CI signal might be varying with a period comparable to the 11- 
year solar cycle. While the evidence for this has weakened, the original claims 
generated renewed interest in neutrino magnetic moment interactions with the 
solar magnetic field. 

The original suggestions by Cisneros and by Okun, Voloshyn, and Vysotsky 
envisioned the rotation 

v e L -> v eR (71) 

producing a right-handed neutrino with sterile interactions in the standard 
model. With the discovery of the MSW mechanism, it was realized that matter 
effects would break the vacuum degeneracy of the v eL and u eR , suppressing the 
spin precession. Lim and MarcianoE3 and AkhmedovL3 pointed out that this 
difficulty was naturally circumvented for 

v e L v m (72) 

as the different matter interactions of the v e and za, can compensate for the 
vacuum v e — mass difference, producing a crossing similar to the usual MSW 
mechanism. Such spin-flavor precession can then occur at full strength due to 
an off-diagonal (in flavor) magnetic moment. 

Quite relevant to this suggestion is the very strong limit on both diagonal 
and off-diagonal magnetic moments is imposed by studies of the red giant 



37 



cooling process of plasmon decay into neutrinos 

7* -n>i&j. (73) 

The result is \fJ,%j\ < 3 ■ 10~ 12 /ib, where hb is an electron Bohr magnetonEl. 
(This can be compared to a simple one-loop estimate c2l of the neutrino mag- 
netic moment of ~ 10~ 1s /ib, taking a typical dark matter value for the neutrino 
mass of a few eV.) If a magnetic moment at the red giant limit is assumed, it 
follows that solar magnetic field strengths of B Q > 10 6 G are needed to produce 
interesting effects. Since the location of the spin-flavor level crossing depends 
on the neutrino energy, such fields have to be extensive to affect an appreciable 
fraction of the neutrino spectrum. It is unclear whether these conditions can 
occur in the sun. 

We will return to a discussion of spin-flavor oscillations and associated 
effects later in these lectures. However, this brief introduction leads naturally 
to the next two topics, neutrino mass possibilities and the evolution of red 
giants. 



3 Dirac and Majorana Neutrinos and Stellar Cooling 

3.1 The Neutrino Mass Matrix 



Consider a general 4n x 4n neutrino mass matrix where n is the number of 
flavors 

/ M L M£\ \ 
Md Ml 


\M* D M R 







Ml Ml 



J 



(74) 



\n/ 

where each entry in this matrix is understood to be a n x n matrix operating in 
flavor space. The entries Mb are the Dirac mass terms, while the Ml and Mr 
are the Majorana terms. The latter break the invariance of the Dirac equation 
under the transformation ipix) — > e %a i\){x) associated with a conserved lepton 
number. Thus it is these terms that govern lepton-number-violating processes 
like double beta decay. 

One can proceed to diagonalize this matrix 



2» 



^v(e) = y U^Viix) with masses nij. 



= 1 



(75) 



The eigenstates are two-component Majorana neutrinos0, yielding the proper 
2 x 2n — An degrees of freedom, where n is the number of flavors. We can 
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recover the Majorana and Dirac limits: 

• If Mr = Ml = 0, the eigenstates of this matrix become pairwise degener- 
ate, allowing the 2n two-component eigenstates to be paired to form n four- 
component Dirac eigenstates. 

• If Mo = 0, the left- and right-handed components decouple, yielding n left- 
handed Majorana eigenstates with standard model interactions. 

There are interesting physical effects associated with these limits. Dirac 
neutrinos can have magnetic dipole, electric dipole (CP and T violating), and 
anapole (P violating) moments, as well as nonzero charge radii. Majorana 
neutrinos can have anapole moments but only transition magnetic and electric 
dipole moments. We mentioned in the previous lecture that transition mag- 
netic moments were quite interesting in the context of MSW effects, as well as 
the use of Mr in the seesaw mechanism. 

In most extended models both Dirac and Majorana mass terms arise, with 
possibilities for interesting physics. We previously described the seesaw mech- 
anism, which explains the lightness of the neutrinos as a effect of Dirac masses 
mixing with a heavy right-handed Majorana mass associated with scales far 
beyond the standard model. Another scenario often discussed in the pseu- 
doDirac limit, where larger Dirac masses are accompanied by small Majorana 
masses. In the limit that the Majorana masses vanish, one recovers the degen- 
eracy that allows one to patch two Majorana neutrinos with opposite CP into 
a four-component Dirac state, as described above. But if the Majorana terms 
remain non- vanishing, two Majorana states remain split by an amount Ml- If 
Ml is much less than Md, one obtains a nearly degenerate pair of states at 
mass Md- 

There are interesting low-energy consequences of these scenarios. Studies 
of the shape of the (3 decay spectrum provides one test of neutrino masses: such 
tests are clearly limited to masses less than the Q- value of the decay. Since the 
different mass eigenstates correspond to distinct final states, the mass eigen- 
states do not interfere. Each neutrino that couples to the electron contributes 
to the decay according to its mixing probability \U e i\ 2 and according to the 
Fermi function F(w e ,m,), where iv e is the electron energy and im the mass 
of the ith neutrino mass eigenstates. Thus massive neutrinos can appear as 
increments to the (3 spectrum, turning on once uj e drops below Q - nii- 

Another interesting process is neutrinoless double (3 decay, which can be 
pictured as (3 decay in which the produced neutrino is reabsorbed on a second 
nucleon, leading to a final state with no neutrinos and two electrons (and thus 
a nuclear charge change of two units). The neutrino mass eigenstates appear 
virtually, and thus clearly interfere. The process is also manifestly lepton- 
numbcr violating, and thus must require Majorana masses. It can be shown 
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that the amplitude, for light neutrino masses, depends on \U e i\ 2 mirjf p , where 
rff p is the relative CP eigenvalue of the ith mass eigenstate. (We are assuming 
CP conservation.) Thus in the Dirac limit, where the mass eigenstates of 
opposite CP become pairwise degenerate, this (3(3 decay mass vanishes. This 
then guarantees that the (3(3 amplitude vanishes in the absence of Majorana 
mass terms. 

3.2 Red Giants and Helium Burning 

We now consider the evolution off the main sequence of a solar-like star, with a 
mass above half a solar mass. As the hydrogen burning in the core progresses 
to the point that no more hydrogen is available, the stellar core consists of 
the ashes from this burning, 4 He. The star then goes through an interesting 
evolution: 

• With no further means of producing energy, the core slowly contracts, thereby 
increasing in temperature as gravity does work on the core. 

• Matter outside the core is still hydrogen rich, and can generate energy 
through hydrogen burning. Thus burning in this shell of material supports 
the outside layers of the star. Note as the core contracts, this matter outside 
the core also is pulled deeper into the gravitational potential. Furthermore, 
the shell H burning continually adds more mass to the core. This means the 
burning in the shell must intensify to generate the additional gas pressure to 
fight gravity. The shell also thickens as this happens, since more hydrogen is 
above the burning temperature. 

• The resulting increasing gas pressure causes the outer envelope of the star 
to expand by a larger factor, up to a factor of 50. The increase in radius 
more than compensates for the increased internal energy generation, so that a 
cooler surface results. The star reddens. Thus this class of star is named a red 
supergiant. 

• This evolution is relatively rapid, perhaps a few hundred million years: the 
dense core requires large energy production. The helium core is supported by 
its degeneracy pressure, and is characterized by densities ~ 10 6 g/cm 3 . This 
stage ends when the core reaches densities and temperatures that allow helium 
burning through the reaction 

a + a + a^ 12 C + 7. (76) 

As this reaction is very temperature dependent (see below), the conditions for 
ignition are very sharply defined. This has the consequence that the core mass 
at the helium flash point is well determined. 

• The onset of helium burning produces a new source of support for the core. 
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The energy release elevates the temperature and the core expands: He burn- 
ing, not electron degeneracy, now supports the core. The burning shell and 
envelope have moved outward, higher in the gravitational potential. Thus shell 
hydrogen burning slows (the shell cools) because less gas pressure is needed 
to satisfy hydrostatic equilibrium. All of this means the evolution of the star 
has now slowed: the red giant moves along the "horizontal branch" , as interior 
temperatures slowly elevate much as in the main sequence. 

The 3a process depends on some rather interesting nuclear physics. The 
first interesting "accident" involves the near degeneracy of the 8 Be ground 
state and two separated as: The 8 Be + ground state is just 92 keV above the 
2a threshold. The measured width of the 8 Be ground state is 2.5 eV, which 
corresponds to a lifetime of 

T m ~ 2.6 • l(T 16 s. (77) 

One can compare this number to the typical time for one a to pass another. 
The red giant core temperature is T 7 ~ 10 -> E ~ 8.6 keV. Thus v/c ~ 0.002. 
So the transit time is 

d 5 f 1 lCT 13 cm 21 

T ~ v ~ 0.002 3 • lO^cm/sec f ( j 

This is more than five orders of magnitude shorter than r m above. Thus when 
a 8 Be nucleus is produced, it lives for a substantial time compared to this naive 
estimate. 

To quantify this, we calculate the flux- averaged cross section assuming 
resonant capture 

{av) •= : { Wr )mT -Y e ~ Er/kT (79) 

where T is the 2a width of the 8 Be ground state. This is the cross section for 
the a + a reaction to form the compound nucleus then decay by a + a. But 
since there is only one channel, this is clearly also the result for producing the 
compound nucleus 8 Be. 

By multiplying the rate/volume for producing 8 Be by the lifetime of 8 Be, 
one gets the number of 8 Be nuclei per unit volume 

N(Be) = 3^(av)T m 
1 + d aa 

N N 1 

= %^r?' 2 e- Er/kT - (80) 
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Notice that the concentration is independent of T. So a small T is not the 
reason we obtain a substantial buildup of 8 Be. This is easily seen: if the width 
is small, then the production rate of 8 Be goes down, but the lifetime of the 
nucleus once it is produced is longer. The two effects cancel to give the same 
8 Be concentration. One sees that the significant 8 Be concentration results from 
two effects: 1) a + a is the only open channel and 2) the resonance energy is 
low enough that some small fraction of the a + a reactions have the requisite 
energy. As E r = 92 keV, E r /kT = 10.67/T 8 (where T 8 is the temperature in 
10 8 K) so that 

N(Be) = A^T~ 3/ V 10 - 67/T8 (0.94 • l(T 33 cm 3 ). (81) 

So plugging in typical values of N a ~ 1.5-10 28 /cm 3 (corresponding to p a ~ 10 5 
g/cm 3 ) and T 8 ~ 1 yields 

N( s Be) 1(1 , , 

iv^T ~ 32 x 10 (82) 

Salpeter suggested that this concentration would then allow a+ 8 Be^ 12 C 
to take place. Hoyle then argued that this reaction would not be fast enough 
to produce significant burning unless it was also resonant. Now the mass of 
8 Be + a, relative to 12 C,is 7.366 MeV, and each nucleus has ,V = 0+. Thus 
s-wave capture would require a + resonance in 12 C at ~ 7.4 MeV. No such 
state was then known, but a search by Cook, Fowler, Lauritsen, and Lauritscn 
revealed a + level at 7.644 MeV, with decay channels 8 Be+a and 7 decay to 
the 2+ 4.433 level in 12 C. The parameters are 

T Q - 8.9eV T 7 - 3.6 • l(T 3 cV. (83) 

The resonant cross section formula gives 

r48 = N8Na{ ^Lf/^ e - Er/kT _ (g4) 

Plugging in our previous expression for iV( 8 Be) yields 

r 48 = 7V 3 T- 3 e- 42 - 9 / Ts (6.3 • 10- 54 cm 6 /s). (85) 
If we denote by iv^ a the decay rate of an a in our plasma, then 
u 5a = 37V 2 T 8 - 3 e- 42 - 9 / Ts (6.3 . 10- 54 cm 6 /sec) 

= ( 1.5.1^ 8 /cm3 ) 2 ( 4 - 3 ' 10 3 /sec)T 8 3 e- 42 -^. (86) 
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Now the energy release per reaction is 7.27 MeV. Thus we can calculate 
the energy produced per gram, e: 

7.27McV 1.5 • 10 23 

£ = ^3a ^ 

3 g 

= (2.5 • ICPerg/g H( 15 ,^ /cm3 )V^. (87) 

We can evaluate this at a temperature of T 8 ~ 1 to find 

e ~ (584ergs/g sec)( - , f ■ (88) 

Typical values found in stellar calculations are in good agreement with this: 
typical red giant energy production is ~ 100 ergs per gram per second. 

To get a feel for the temperature sensitivity of this process, we can do a 
Taylor series expansion, finding 

e(T) ~ (^r o Nl (89) 

- 1 o 

This steep temperature dependence is the reason the He flash is delicately de- 
pendent on conditions in the core. 



3.3 Neutrino Magnetic Moments and He Ignition 

Prior to the helium flash, the degenerate He core radiates energy largely by 
neutrino pair emission. The process is the decay of a plasmon — which one 
can think of as a photon "dressed" by electron-hole excitations — thereby 
acquiring an effective mass of about 10 keV. The photon couples to a neutrino 
pair through a electron particle-hole pair that then decays into aZ„-> vv. 

If this cooling is somehow enhanced, the degenerate helium core would 
be kept cooler, and would not ignite at the normal time. Instead it would 
continue to grow until it overcame the enhanced cooling to reach, once again, 
the ignition temperature. 

One possible mechanism for enhanced cooling is a neutrino magnetic mo- 
ment. Then the plasmon could directly couple to a neutrino pair. The strength 
of this coupling would depend on the size of the magnetic moment. 

A delay in the time of He ignition has several observable consequences, 
including changing the ratio of red giant to horizontal branch stars. Thus, using 
the standard theory of red giant evolution, investigators have attempted to 
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determine what size of magnetic momept would produce unacceptable changes 
in the astronomy. The result is a limit cB on the neutrino magnetic moment of 

Hij < 3 • 10~ 12 clcctron Bohr magnetons (90) 

as was mentioned earlier. This limit is more than two orders of magnitude 
more stringent than that from direct laboratory tests. 

This example is just one of a number of such constraints that can be 
extracted from similar stellar cooling arguments. The arguments above, for 
example, can be repeated for neutrino electric dipole moments. More inter- 
esting, it can be repeated for axion emission from red giants. Axions, the 
pseudoGolstone bosons associate with the solution of the strong CP problem 
suggested by Peccei and Quinn, are very light and can be produced radiatively 
within the red giant by the Compton process, by the Primakoff process off 
nuclei, or by emission from low-lying nuclear levels, such at from the 14 keV 
transition in 57 Fe. The net result is that axions of mass above a few eV are 
excluded; if axions have a direct coupling to electrons, so that the Compton 
process off electrons operates, the constraint is considerably tighter. 

A similar argument can be formulated for supernova cooling. During 
SN1987A the neutrino burst detected by 1MB and by Kamiokande was consis- 
tent with cooling on a timescale of about 4 seconds. Thus any process cooling 
the star more efficiently than neutrino emission would have shortened this time, 
while also reducing the flux in neutrinos. Large Dirac neutrino masses allow 
trapped neutrinos to scatter into sterile right-handed states. Right-handed 
neutrinos, lacking standard model interactions, would then escape the star 
(provided they do not scatter back into interacting left-handed states). Un- 
fortunately the upper bounds imposed on the neutrino mass are quite model 
dependent, ranging over (1-25) keV. 

The supernova cooling argument can also be repeated for axions. The 
window of sensitive runs from 1 eV (above this mass they are more strongly 
coupled than neutrinos, and thus cannot compete with neutrino cooling) to 
about 0.01 eV (below this mass they are too weakly coupled to be produced 
on the timescale of supernova cooling). It is interesting that the supernova 
and red giant cooling limit on axions nearly meet: a small window may still 
exist around a few eV if the axion has no coupling to electrons. 

4 Atmospheric Neutrinos 

Atmospheric neutrinos arise from the decay of secondary pions, kaons, and 
muons produced by the collisions of primary cosmic rays with the oxygen and 
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nitrogen nuclei in the upper atmosphere. For energies less than 1 GeV all the 
secondaries decay : 

->e ± + i/ e (F e )+F A1 (^ I ). (91) 

Consequently one expects the ratio 

r = (t/ e + V e )/(v lt +V IA ) (92) 

to ba-approximately 0.5 in this energy range. Detailed Monte Carlo calcula- 
tions c3, including the effects of muon polarization, give r ~ 0.45. Since one 
is evaluating a ratio of similarly calculated processes, systematic errors are 
significantly reduced. Different groups estimating this ratio, even though they 
start with neutrino fluxes which can differ in magnitude by up to 25%, all agree 
within a few percenter. As the shower energy increases more muons survive 
due to time dilation. Hence one expects the ratio r to decrease as the energy 
increases. The ratio (observed to predicted) of ratios 

R = (^/"eWa (Q3) 
(^/J / e)MonteCarlo 

was determined in several experiments. There seems to be a persistent dis- 
crepancy between theory and experiment. The most recent measurement of 
this ratio at SupcrKamiokandeE-S gives 

R = 0.63 ± 0.03(stat) ± 0.05(syst) (94) 

for sub-GeV events which were fully contained in the detector and 

R = 0.65 ± 0.05(stat) ± 0.08(syst) (95) 

for fully- and partially-contained multi-GcV events. As the agreement between 
data and theoretical expectation would give R — 1, neutrino oscillations are in- 
voked to explain the discrepancy. Another evidence for this explanation comes 
from measuring R as a function of the zenith angle, 0, between the vertical 
and neutrino direction. A down-going neutrino (0 ~ 0°) travels through the 
atmosphere above the detector (a distance of about 20 km), whereas an up- 
going neutrino (0 ~ 180°) has traveled through the entire Earth (a distance of 
about 13000 km). Hence a measurement of number of neutrinos as a function 
of the zenith angle yields information about their numbers as a function of the 
distance traveled. 
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Figure 13: The ratio of fully contained events measured at SuperKamiokande versus recon- 
structed L/Ev. the dashed lines show the expected shape for Vp <-> v T oscillations with 
Sm 2 = 2.2 ■ 1CT 3 eV 2 and sin 2 29 = 1. 



The SuperKamiokande collaboration measured the zenith angle depen- 
nce of not only R, but of the electron and muon neutrino fluxes separately 
. This information is shown in Figure 13, where the data is plotted as a 
function of the reconstructed L/E v instead of the zenith angle. The data ex- 
hibit a zenith-angle (distance) dependent deficit of muon neutrinos, but not of 
electron neutrinos. This behavior is consistent with <-> v T oscillations. This 
interpretation is-consistent with -the recent measurements of up-going muons 
at Kamiokandec3 and MACRO c3 detectors. These muons are produced by 
the very high energy up-going muon neutrinos in the rocks surrounding the 
detectors. Both of these experiments also observe a deficit. 

5 Supernovae, Supernova Neutrinos, and Nucleosynthesis 

Consider a massive star, in excess of 10 solar masses, burning the hydrogen in 
its core under the conditions of hydrostatic equilibrium. When the hydrogen 
is exhausted, the core contracts until the density and temperature are reached 
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where 3a -^ 12 G can take place. The He is then burned to exhaustion. This 
pattern (fuel exhaustion, contraction, and ignition of the ashes of the previous 
burning cycle) repeats several times, leading finally to the explosive burning 
of 28 Si to Fe. For a heavy star, the evolution is rapid: the star has to work 
harder to maintain itself against its own gravity, and therefore consumes its 
fuel faster. A 25 solar mass star would go through all of these cycles in about 
7 My, with the final explosive Si burning stage taking a few days. The result 
is an "onion skin" structure of the pre-collapse star in which the star's history 
can be read by looking at the surface inward: there are concentric shells of H, 
4 He, 12 C, 16 and 20 Ne, 28 Si, and 56 Fe at the center. 

5.1 The Explosion Mechanism^^ 

The source of energy for this evolution is nuclear binding energy. A plot of the 
nuclear binding energy 5 as a function of nuclear mass shows that the minimum 
is achieved at Fe. In a scale where the 12 C mass is picked as zero: 
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Thus once the Si burns to produce Fe, there is no further source of nuclear 
energy adequate to support the star. So as the last remnants of nuclear burning 
take place, the core is largely supported by degeneracy pressure, with the 
energy generation rate in the core being less than the stellar luminosity. The 
core density is about 2 xlO 9 g/cc and the temperature is kT ~ 0.5 MeV. 

Thus the collapse that begins with the end of Si burning is not halted 
by a new burning stage, but continues. As gravity does work on the matter, 
the collapse leads to a rapid heating and compression of the matter. As the 
nucleons in Fe are bound by about 8 MeV, sufficient heating can release as and 
a few nucleons. At the same time, the electron chemical potential is increasing. 
This makes electron capture on nuclei and any free protons favorable, 

e~ + p — > v e + n. (96) 

Note that the chemical equilibrium condition is 

Me + = + (97) 

Thus the fact that neutrinos are not trapped plus the rise in the electron Fermi 
surface as the density increases, lead to increased neutronization of the matter. 
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The escaping neutrinos carry off energy and lepton number. Both the electron 
capture and the nuclear excitation and disassociation take energy out of the 
electron gas, which is the star's only source of support. This means that the 
collapse is very rapid. Numerical simulations find that the iron core of the star 
(~ 1.2-1.5 solar masses) collapses at about 0.6 of the free fall velocity. 

In the early stages of the infall the v e s readily escape. But neutrinos 
are trapped when a density of ~ 10 12 g/cm 3 is reached. At this point the 
neutrinos begin to scatter off the matter through both charged current and 
coherent neutral current processes. The neutral current neutrino scattering 
off nuclei is particularly important, as the scattering cross section is off the 
total nuclear weak charge, which is approximately the neutron number. This 
process transfers very little energy because the mass energy of the nucleus is 
so much greater than the typical energy of the neutrinos. But momentum is 
exchanged. Thus the neutrino "random walks" out of the star. When the 
neutrino moan free path becomes sufficiently short, the "trapping time" of the 
neutrino begins to exceed the time scale for the collapse to be completed. This 
occurs at a density of about 10 12 g/cm 3 , or somewhat less than 1% of nuclear 
density. After this point, the energy released by further gravitational collapse 
and the star's remaining lepton number are trapped within the star. 

If we take a neutron star of 1.4 solar masses and a radius of 10 km, an 
estimate of its binding energy is 

CM 2 

— — ~ 2.5 x 10 53 ergs. (98) 
2R 

Thus this is roughly the trapped energy that will later be radiated in neutrinos. 

The trapped lepton fraction Yj, is a crucial parameter in the explosion 
physics: a higher trapped Yl leads to a larger homologous core, a stronger 
shock wave, and easier passage of the shock wave through the outer core, as 
will be discussed below. Most of the lepton number loss of an infalling mass 
element occurs as it passes through a narrow range of densities just before 
trapping. The reasons for this are relatively simple: on dimensional grounds 
weak rates in a plasma go as T 5 , where T is the temperature. Thus the electron 
capture rapidly turns on as matter falls toward the trapping radius, and lepton 
number loss is maximal just prior to trapping. Inelastic neutrino reactions have 
an important effect on these losses, as the coherent trapping cross section goes 
as E"l and is thus least effective for the lowest energy neutrinos. As these 
neutrinos escape, inelastic reactions repopulate the low energy states, allowing 
the neutrino emission to continue. 

The velocity of sound in matter rises with increasing density. The inner 
homologous core, with a mass Mhc ^0.6 — 0.9 solar masses, is that part of 
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the iron core where the sound velocity exceeds the infall velocity. This allows 
any pressure variations that may develop in the homologous core during infall 
to even out before the collapse is completed. As a result, the homologous core 
collapses as a unit, retaining its density profile. That is, if nothing were to 
happen to prevent it, the homologous core would collapse to a point. 

The collapse of the homologous core continues until nuclear densities are 
reached. As nuclear matter is rather incompressible (~ 200 MeV/f 3 ), the 
nuclear equation of state is effective in halting the collapse: maximum densities 
of 3-4 times nuclear are reached, e.g., perhaps 6 • 10 14 g/cm 3 . The innermost 
shell of matter reaches this supernuclear density first, rebounds, sending a 
pressure wave out through the homologous core. This wave travels faster than 
the infalling matter, as the homologous core is characterized by a sound speed 
in excess of the infall speed. Subsequent shells follow. The resulting series of 
pressure waves collect near the sonic point (the edge of the homologous core). 
As this point reaches nuclear density and comes to rest, a shock wave breaks 
out and begins its traversal of the outer core. 

Initially the shock wave may carry an order of magnitude more energy 
than is needed to eject the mantle of the star (less than 10 51 ergs). But 
as the shock wave travels through the outer iron core, it heats and melts 
the iron that crosses the shock front, at a loss of ~ 8 MeV/nucleon. The 
enhanced electron capture that occurs off the free protons left in the wake of 
the shock, coupled with the sudden reduction of the neutrino opacity of the 
matter (recall a co herent ~ N 2 ), greatly accelerates neutrino emission. This is 
another energy loss]] The summed losses from shock wave heating and neutrino 
emission are comparable to the initial energy carried by the shock wave. Thus 
most numerical models fail to produce a successful "prompt" hydrodynamic 
explosion. 

Most of the attention in the past decade focused on two explosion scenarios. 
In the prompt mechanism described above, the shock wave is sufficiently strong 
to survive the passage of the outer iron core with enough energy to blow off 
the mantle of the star. The most favorable results were achieved with smaller 
stars (less than 15 solar masses) where there is less overlying iron, and with soft 
equations of state, which produce a more compact neutron star and thus lead 
to more energy release. In part because of the lepton number loss problems 
discussed earlier, now it is widely believed that this mechanism fails for all but 
unrealistically soft nuclear equations of state. 

b Many numerical models predict a strong "breakout" burst of v e s in the few milliseconds 
required for the shock wave to travel from the edge of the homologous core to the neutri- 
nosphere at p ~ 10 12 g/cm 3 and r ~ 50 km. The neutrinosphcrc is the term from the 
neutrino trapping radius, or surface of last scattering. 
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The delayed mechanism begins with a failed hydrodynamic explosion; after 
about 0.01 seconds the shock wave stalls at a radius of 200-300 km. It exists in 
a sort of equilibrium, gaining energy from matter falling across the shock front, 
but loosing energy to the heating of that material. However, after perhaps 0.5 
seconds, the shock wave is revived due to neutrino heating of the nucleon 
"soup" left in the wake of the shock. This heating comes primarily from 
charged current reactions off the nucleons in that nucleon gas; quasi-elastic 
scattering also contributes. This high entropy radiation-dominated gas may 
reach two McV in temperature. The pressure exerted by this gas helps to push 
the shock outward. It is important to note that there are limits to how effective 
this neutrino energy transfer can be: if matter is too far from the core, the 
coupling to neutrinos is too weak to deposit significant energy. If too close, the 
matter may be at a temperature (or soon reach a temperature) where neutrino 
emission cools the matter as fast or faster than neutrino absorption heats it. 
The term "gain radius" is used to describe the region where useful heating is 
done. 

This subject is still controversial and unclear. The problem is numerically 
challenging, forcing modelers to handle the difficult hydrodynamics of a shock 
wave; the complications of the nuclear equation of state at densities not yet 
accessible to experiment; modeling in two or three dimensions; handling the 
slow diffusion of neutrinos; etc. Not all of these aspects can be handled rea- 
sonably at the same time, even with existing supercomputers. Thus there is 
considerable disagreement about whether we have any supernova model that 
succeeds in ejecting the mantle. 

However the explosion proceeds, there is agreement that 99% of the 3 -10 53 
ergs released in the collapse is radiated in neutrinos of all flavors. The time 
scale over which the trapped neutrinos leak out of the protoneutron star is 
about three seconds. Through most of their migration out of the protoneutron 
star, the neutrinos are in flavor equilibrium 

e.g., v c + D c <-> + v^. (99) 

As a result, there is an approximate cquipartition of energy among the neutrino 
flavors. After weak decoupling, the v e s and i7 e s remain in equilibrium with the 
matter for a longer period than their heavy-flavor counterparts, due to the 
larger cross sections for scattering off electrons and because of the charge- 
current reactions 

v e + n <-> p + e~ 

v e +p <-> n + e + . (100) 
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Thus the heavy flavor neutrinos decouple from deeper within the star, where 
temperatures are higher. Typical calculations yield 

T Ufl ~ T Vt ~ 8MeV T„ e ~ 3.5MeV TV e - 4.5MeV. (101) 

The difference between the v e and v e temperatures is a result of the neutron 
richness of the matter, which enhances the rate for charge-current reactions of 
the v e s, thereby keeping them coupled to the matter somewhat longer. 

This temperature hierarchy is crucially important to nucleosynthesis and 
also to possible neutrino oscillation scenarios. The three- flavor MSW level- 
crossing diagram is shown in Fig. 14. One very popular scenario attributes 
the solar neutrino problem to <-* v e transmutation; this means that a sec- 
ond crossing with a v T could occur at higher density. It turns out plausible 
seesaw mass patterns suggest a v r mass on the order of a few eV, which would 
be interesting cosmologically. The second crossing would then occur outside 
the neutrino sphere, that is, after the neutrinos have decoupled and have fixed 
spectra with the temperatures given above. Thus a v e «-> v T oscillation would 
produce a distinctive T ~ 8 MeV spectrum of u e s. This has dramatic conse- 
quences for terrestrial detection and for nucleosynthesis in the supernova. 

5.2 The Neutrino Process^ 

Core-collapse supernovae are one of the major engines driving galactic chem- 
ical evolution, producing and ejecting the metals that enrich our galaxy. The 
discussion of the previous section described the hydrostatic evolution of a pre- 
supernova star in which large quantities of the most abundant metals (C, O, 
Ne,...) are synthesized and later ejected during the explosion. During the pas- 
sage of the shock wave through the star's mantle, temperature of ~ (1 — 3) • 10 9 K 
and are reached in the silicon, oxygen, and neon shells. This shock wave heat- 
ing induces (7, a) «-* (a, 7) and related reactions that generate a mass flow 
toward highly bound nuclei, resulting in the synthesis of iron peak elements 
as well as less abundant odd-A species. Rapid neutron-induced reactions are 
thought to take place in the high-entropy atmosphere just above the mass cut, 
producing about half of the heavy elements above A ^ 80. This is the sub- 
ject of the Sec. 4.3. Finally, the i/-process described below is responsible for 
the synthesis of rare species such as 11 B and 19 F. This process involves the 
response of nuclei at momentum transfers where the allowed approximation 
is no longer valid. Thus we will use the ^-process in this section to illustrate 
some of the relevant nuclear physics. 

One of the problems - still controversial - that may be connected with the 
neutrino process is the origin of the light elements Be, B and Li, elements which 
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Figure 14: Three-flavor neutrino level-crossing diagram. One popular scenario associates the 
solar neutrino problem with v a <-» oscillations and predicts a cosmologically interested 
massive v T with v e <-> v T oscillations near the supernova neutrinosphcrc. 



are not produced in sufficient amounts in the big bang or in any of the stellar 
mechanisms we have discussed. The traditional explanation has been cosmic 
ray spallation interactions with C, O, and N in the interstellar medium. In this 
picture, cosmic ray protons collide with C at relatively high energy, knocking 
the nucleus apart. So in the debris one can find nuclei like 10 B, n B, and 7 Li. 

But there are some problems with this picture. First of all, this is an 
example of a secondary mechanism: the interstellar medium must be enriched 
in the C, O, and N to provide the targets for these reactions. Thus cosmic ray 
spallation must become more effective as the galaxy ages. The abundance of 
boron, for example, would tend to grow quadratically with metallicity, since 
the rate of production goes linearly with metallicity. But observations, espe- 
cially recent measurements with the HST, find a linear growthca in the boron 
abundance. 

A second problem is that the spectrum of cosmic ray protons peaks near 
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1 GeV, leading to roughly comparable production of the two isotopes 10 B 
and n B. That is, while it takes more energy to knock two nucleons out of 
carbon than one, this difference is not significant compared to typical cosmic 
ray energies. More careful studies lead to the expectation that the abundance 
ratio of n B to 10 B might be ~ 2. In nature, it is greater than 4. 

Fans of cosmic ray spallation have offered solutions to these problems, 
e.g., similar reactions occurring in the atmospheres of nebulae involving lower 
energy cosmic rays. As this suggestion was originally stimulated by the ob- 
servation of nuclear 7 rays from Orion, now retracted, some of the motivation 
for this scenario has evaporated. Here we focus on an alternative explanation, 
synthesis via neutrino spallation. 

Previously we described the allowed Gamow- Teller (spin-flip) and Fermi 
weak interaction operators. These are the appropriate operators when one 
probes the nucleus at a wavelength that is, at a size scale - where the nucleus 
responds like an elementary particle. We can then characterize its response by 
its macroscopic quantum numbers, the spin and charge. On the other hand, 
the nucleus is a composite object and, therefore, if it is probed at shorter length 
scales, all kinds of interesting radial excitations will result, analogous to the 
vibrations of a drumhead. For a reaction like neutrino scattering off a nucleus, 
the full operator involves the additional factor 

e llf ~l + ik-r (102) 

where the expression on the right is valid if the magnitude of k is not too large. 
Thus the full charge operator includes a "first forbidden" term 

A 

J2nr 3 (i) (103) 

i=i 

and similarly for the spin operator 

A 

J]Pi®5(t)]j=o,i,2is(i). (104) 

These operators generate collective radial excitations, leading to the so-called 
"giant resonance" excitations in nuclei. The giant resonances are typically at 
an excitation energy of 20-25 MeV in light nuclei. One important property is 
that these operators satisfy a sum rule (Thomas-Reichc-Kuhn) of the form 

EK/lE r « T 3(0K>l 2 ~™ \ (105) 

/ »=i 
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where the sum extends over a complete set of final nuclear states. These 
first-forbidden operators tend to dominate the cross sections for scattering the 
high energy supernova neutrinos (v^s and z^s), with E v ~ 25 MeV, off light 
nuclei. From the sum rule above, it follows that nuclear cross sections per 
target nucleon are roughly constant. 

The Ef giant dipole mode described above is depicted qualitatively in Fig. 
15a. This description, which corresponds to an early model of the giant res- 
onance response by Goldhaber and Teller, involves the harmonic oscillation 
of the proton and neutron fluids against one another. The restoring force for 
small displacements would be linear in the displacement and dependent on the 
nuclear symmetry energy. There is a natural extension of this model to weak 
interactions, where axial excitations occur. For example, one can envision a 
mode similar to that of Fig. 15a where the spin-up neutrons and spin-down pro- 
tons oscillate against spin-down neutrons and spin-up protons, the spin-isospin 
mode of Fig. 15b. This mode is one that arises in a simple SU(4) extension of 
the Goldhaber- Teller model, derived by assuming that the nuclear force is spin 
and isospin independent, at the same excitation energy as the El mode. In 
full, the Goldhaber- Teller model predicts a degenerate 15-dimcnsional super- 
multiplet of giant resonances, each obeying sum rules analogous to the TRK 
sum rule. While more sophisticated descriptions of the giant resonance region 
are available, of course, this crude picture is qualitatively accurate. 




Figure 15: Schematic illustration of a) the El giant dipole mode familiar from electromag- 
netic interactions and b) a spin-isospin giant dipole mode associated with the first-forbidden 

weak axial response. 

This nuclear physics is important to the ^-process. The simplest example 
of i^-process nucleosynthesis involves the Ne shell in a supernova. Because of 
the first-forbidden contributions, the cross section for inelastic neutrino scat- 
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tcring to the giant resonances in Ne is <~ 3 • 10 41 cm 2 /flavor for the more 
energetic heavy-flavor neutrinos. This reaction 

v + A -> v' + A* (106) 

transfers an energy typical of giant resonances, ~ 20 McV. A supernova releases 
about 3 x 10 53 ergs in neutrinos, which converts to about 4 x 10 57 heavy flavor 
neutrinos. The Ne shell in a 20 M© star has at a radius ~ 20,000 km. Thus 
the neutrino fluence through the Ne shell is 

~ 10 38 /cm 2 . (107) 



r 4tt(20, 000km) 2 

Thus folding the fluence and cross section, one concludes that approximately 
l/300th of the Ne nuclei interact. 

This is quite interesting since the astrophysical origin of 19 F had not been 
understood. The only stable isotope of fluorine, 19 F has an abundance 

l»F 1 

(108) 



20 Ne 3100' 

This leads to the conclusion that the fluorine found in a tube of toothpaste 
was created by neutral current neutrino reactions deep inside some ancient 
supernova. 

The calculation of the final 19 F/ 20 Ne ratio is more complicated than the 
simple 1/300 ratio given above: 

• When Ne is excited by ~ 20 McV through inelastic neutrino scattering, it 
breaks up in two ways 

20 Ne(^, z/) 20 Ne* Ne + n ^ 19 F + e+ + ^ c + n 
20 Nc(^y) 2 °Ne* ^ 19 F + p (109) 

with the first reaction occurring half as frequently as the second. As both 
channels lead to 19 F, we have correctly estimated the instantaneous abundance 
ratio in the Ne shell of 

i9 F i 

*Ne ~ 300- (110) 

• We must also address the issue of whether the produced 19 F survives. In the 
first 10 -8 sec the co-produced neutrons in the first reaction react via 

15 0(n,p) 15 N 19 Ne(n,a) 16 20 Ne(n, 7 ) 21 Ne 19 Ne(n,p) 19 F (111) 
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with the result that about 70% of the 19 F produced via spallation of neutrons 
is then immediate destroyed, primarily by the (n, a) reaction above. In the 
next 10 -6 sec the co-produced protons are also processed 

15 N(p,a) 12 C 19 F(p,a) 16 23 Na(p, a) 20 Ne (112) 

with the latter two reactions competing as the primary proton poisons. This 
makes an important prediction: stars with high Na abundances should make 
more F, as the 23 Na acts as a proton poison to preserve the produced F. 
• Finally, there is one other destruction mechanism, the heating associated 
with the passage of the shock wave. It turns out the the F produced prior to 
shock wave passage can survive if it is in the outside half of the Ne shell. The 
reaction 

19 F( 7 ,a) 15 N (113) 

destroys F for peak explosion temperatures exceeding 1.7-10 9 K. Such a temper- 
ature is produced at the inner edge of the Ne shell by the shock wave heating, 
but not at the outer edge. 

If all of this physics in handled is a careful network code that includes 
the shock wave heating and F production both before and after shock wave 
passage, the following are the results: 

[ 19 F/ 20 Ne]/[ 19 F/ 20 Ne] e T hoavy y (MeV) 
0.14 4 
0.6 6 
1.2 8 
1.1 10 
1.1 12 



where the abundance ratio in the first column has been normalized to the solar 
value. One sees that the attribution of F to the neutrino process argues that 
the heavy flavor v temperature must be greater than 6 MeV, a result theory 
favors. One also sees that F cannot be overproduced by this mechanism: 
although the instantaneous production of F continues to grow rapidly with 
the neutrino temperature, too much F results in its destruction through the 
(p, a) reaction, given a solar abundance of the competing proton poison 23 Na. 
Indeed, this illustrates an odd quirk: although in most cases the neutrino 
process is a primary mechanism, one needs 23 Na present to produce significant 
F. Thus in this case the neutrino process is a secondary mechanism. 

While there are other significant neutrino process products ( 7 Li, 138 La, 
180 Ta, 15 N ...), the most important product is n B, produced by spallation off 
carbon. A calculation by Timmes et alu3 found that the combination of the 
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neutrino process, cosmic ray spallation and big-bang nucleosynthesis together 
can explain the evolution of the light elements. The neutrino process, which 
produces a great deal of 11 B but relatively little 10 B, combines with the cos- 
mic ray spallation mechanism to yield the observed isotope ratio. Again, one 
prediction of this picture is that early stars should be n B rich, as the neutrino 
process is primary and operates early in our galaxy's history; the cosmic ray 
production of 10 B is more recent. There is hope that HST studies will soon be 
able to discriminate between 10 B and n B: as yet this has not been done. 

5.3 The r-process 

Beyond the iron peak nuclear Coulomb barriers become so high that charged 
particle reactions become ineffective, leaving neutron capture as the mechanism 
responsible for producing the heaviest nuclei. If the neutron abundance is 
modest, this capture occurs in such a way that each newly synthesized nucleus 
has the opportunity to (3 decay, if it is energetically favorable to do so. Thus 
weak equilibrium is maintained within the nucleus, so that synthesis is along 
the path of stable nuclei. This is called the s- or slow-process. However a plot 
of the s-process in the (N,Z) plane reveals that this path misses many stable, 
neutron-rich nuclei that are known to exist in nature. This suggests that 
another mechanism is at work, too. Furthermore, the abundance peaks found 
in nature near masses A ~ 130 and A ~ 190, which mark the closed neutron 
shells where neutron capture rates and (3 decay rates are slower, each split into 
two sub-peaks. One set of sub-peaks corresponds to the closed-neutron-shell 
numbers N ~ 82 and N ~ 126, and is clearly associated with the s-process. The 
other set is shifted to smaller N, <~ 76 and ~ 116, respectively, and is suggestive 
of a much more explosive neutron capture environment where neutron capture 
can be rapid. 

This second process is the r- or rapid-process, characterized by: 

• The neutron capture is fast compared to (3 decay rates. 

• The equilibrium maintained within a nucleus is established by (n, 7) <-» (7, n): 
neutron capture fills up the available bound levels in the nucleus until this 
equilibrium sets in. The new Fermi level depends on the temperature and the 
relative n/7 abundance. 

• The nucleosynthesis rate is thus controlled by the (3 decay rate: each (3~ 
capture converting n — > p opens up a hole in the neutron Fermi sea, allowing 
another neutron to be captured. 

• The nucleosynthesis path is along exotic, neutron-rich nuclei that would be 
highly unstable under normal laboratory conditions. 

• As the nucleosynthesis rate is controlled by the j3 decay, mass will build up 
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at nuclei where the (3 decay rates are slow. It follows, if the neutron flux is 
reasonable steady over time so that equilibrated mass flow is reached, that the 
resulting abundances should be inversely proportional to these (3 decay rates. 

Let's first explore the (n, 7) ^> (7, n) equilibrium condition, which requires 
that the rate for (71,7) balances that for (7,71) for an average nucleus. So 
consider the formation cross section 

A + n^ (114) 

This is an exothermic reaction, as the neutron drops into the nuclear well. Our 
averaged cross section, assuming a resonant reaction (the level density is high 
in heavy nuclei) is 

where E <~ is the resonance energy, and the Ts are the indicated partial and 
total widths. Thus the rate per unit volume is 

( 2tt \ 3/2 r„r 7 



^-N n N A ^—) (116) 



where N n and Na are the neutron and nuclear number densities and /i the 
reduced mass. This has to be compared to the (7, n) rate. 

The (7, n) reaction requires the photon number density in the gas. This 
is given by the Bose-Einstein distribution 

at/ \ 87T e 2 de 
The high-energy tail of the normalized distribution can thus be written 



' 2 --*/ kT de (118) 



iV 7 7r 2 



e e 



where in the last expression we have set % = c = 1. 

Now we need the resonant cross section in the (7, n) direction. For photons 
the wave number is proportional to the energy, so 

CT(7 ' n) = ? (e - E r f + (r/2) 2 ' (1 } 

As the velocity is c =1, 

<">'i^f (120) 
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We evaluate this in the usual way for a sharp resonance, remembering that the 
energy integral over just the denominator above (the sharply varying part) is 
27r/r 

~ ElEUe^"/^-. (121) 



N 7 



So that the rate becomes 



r {l>n) ~2N A+ ^e- E ^ T . (122) 
Equating the (n, 7) and (7, n) rates and taking Na ~ Na-i then yields 

N »~jk{^T e - B ' itT <i23 » 

where the fis and cs have been properly inserted to give the right dimensions. 
Now E r is essentially the binding energy. So plugging in the conditions N n ~ 
3 x 10 23 /cm 3 and T 9 ~ 1, we find that the binding energy is ~ 2.4 MeV. Thus 
neutrons are bound by about 30 times kT, a value that is still small compared 
to a typical binding of 8 MeV for a normal nucleus. (In this calculation the 
neutron reduced mass is calculated by assuming a nuclear target with A=150.) 

The above calculation fails to count spin states for the photons and nuclei 
and is thus not quite correct. But it makes the essential point: the r-process 
involves very exotic species largely unstudied in any terrestrial laboratory. It 
is good to bear this in mind, as in the following section we will discuss the 
responses of such nuclei to neutrinos. Such responses thus depend on the 
ability of theory to extrapolate responses from known nuclei to those quite 
unfamiliar. 

The path of the r-process is along neutron-rich nuclei, where the neutron 
Fermi sea is just <~ (2-3) MeV away from the neutron drip line (where no 
more bound neutron levels exist). After the r-process finishes (the neutron 
exposure ends) the nuclei decay back to the valley of stability by (3 decay. This 
can involve some neutron spallation (/3-delayed neutrons) that shift the mass 
number A to a lower value. But it certainly involves conversion of neutrons 
into protons, and that shifts the r-process peaks at N ~ 82 and 126 to a lower 
N, off course. This effect is clearly seen in the abundance distribution: the 
r-process peaks are shifted to lower N relative to the s-process peaks. This 
is the origin of the second set of "sub-peaks" mentioned at the start of the 
section. 

It is believed that the r-process can proceed to very heavy nuclei (A ~ 
270) where it is finally ended by /3-delayed and n-induced fission, which feeds 
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matter back into the process at an A ~ A max /2. Thus there may be important 
cycling effects in the upper half of the r-process distribution. 

What is the site(s) of the r-process? This has been debated many years 
and still remains a controversial subject: 

• The r-process requires exceptionally explosive conditions 

p(n) - 10 20 cm" 3 T - 10 9 K t - Is. 

• Both primary and secondary sites proposed. Primary sites are those not 
requiring preexisting metals. Secondary sites are those where the neutron cap- 
ture occurs on preexisting s-process seeds. 

• Suggested primary sites include the the neutronized atmosphere above the 
proto-neutron star in a Type II supernova, neutron-rich jets produced in su- 
pernova explosions or in neutron star mergers, inhomogeneous big bangs, etc. 

• Secondary sites, where p(n) can be lower for successful synthesis, include the 
He and C zones in Type II supernovae, the red giant He flash, etc. 

The balance of evidence favors a primary site, so one requiring no pre- 
cnrichmcnt of heavy s-process metals. Among the evidence: 

1) HST studies of very- metal-poor halo stars: The most important evidence 
are the recent HST measurements of Cowan, Sneden et al. E2l of very metal- 
poor stars ([Fe/H] ~ -1.7 to -3.12) where an r-process distribution very much 
like that of our sun has been seen for Z > 56. Furthermore, in these stars 
the iron content is variable. This suggests that the "time resolution" inher- 
ent in these old stars is short compared to galactic mixing times (otherwise 
Fe would be more constant). The conclusion is that the r-process material 
in these stars is most likely from one or a few local supernovae. The fact 
that the distributions match the solar r-process (at least above charge 56) 
strongly suggests that there is some kind of unique site for the r-process: the 
solar r-process distribution did not come from averaging over many different 
kinds of r-process events. Clearly the fact that these old stars are enriched 
in r-process metals also strongly argues for a primary process: the r-process 
works quite well in an environment where there are few initial s-process metals. 

2) There are also fairly good theoretical arguments that a primary r-process oc- 
curring in a core-collapse supernova might be viableEd. First, galactic chemical 
evolution studies indicate that the growth of r-process elements in the galaxy 
is consistent with low-mass Type II supernovae in rate and distribution. More 
convincing is the fact that modelers have shown that the conditions needed for 
an r-process (very high neutron densities, temperatures of 1-3 billion degrees) 
might be realized in a supernova. The site is the last material expelled from 
the supernova, the matter just above the mass cut. When this material is 
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blown off the star initially, it is a very hot neutron-rich, radiation-dominated 
gas containing neutrons and protons, but an excess of the neutrons. As it 
expands off the star and cools, the material first goes through a freeze-out to 
a particles, a step that essentially locks up all the protons in this way. Then 
the as interact through reactions like 

a + a + a -> 12 C 
a + a + n — > 9 Be 

to start forming heavier nuclei. Note, unlike the big bang, that the density 
is high enough to allow such three-body interactions to bridge the mass gaps 
at A = 5,8. The a capture continues up to heavy nuclei, to A ~ 80, in the 
network calculations. The result is a small number of "seed" nuclei, a large 
number of as, and excess neutrons. These neutrons preferentially capture on 
the heavy seeds to produce an r-process. Of course, what is necessary is to 
have ~ 100 excess neutrons per seed in order to successfully synthesize heavy 
mass nuclei. Some of the modelers find conditions where this almost happens. 

There are some very nice aspects of this site: the amount of matter ejected 
is about 10~ 5 — 10~ 6 solar masses, which is just about what is needed over the 
lifetime of the galaxy to give the integrated r-process metals we see, taking a 
reasonable supernova rate. But there are also a few problems, especially the 
fact that with calculated entropies in the nucleon soup above the proto-neutron 
star, neutron fractions appear to be too low to produce a successful A ~ 190 
peak. There is some interesting recent work invoking neutrino oscillations 
to cure this problem: charge current reactions on free protons and neutrons 
determine the n/p ratio in the gas. Then, for example, an oscillation of the 
type v e — > ^sterile can alter this ratio, aiit would turn off the v e s that destroy 
neutrons by charged-current reactions^. Unfortunately, a full discussion of 
such possibilities would take us too far afield today. 

The nuclear physics of the r-process tells us that the synthesis occurs when 
the nucleon soup is in the temperature range of (3-1) T0 9 K, which, in the hot 
bubble r-process described above, corresponds to a freeze-out radius of (600- 
100) km and a time ~ 10 seconds after core collapse. The neutrino fluence after 
freeze-out (when the temperature has dropped below 10 9 K and the r-process 
stops) is then ~ (0.045-0.015) -10 51 ergs/(100km). Thus, after completion of 
the r-process, the newly synthesized material experiences an intense flux of 
neutrinos. This brings up the question of whether the neutrino flux could have 
any effect on the r-process. 
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5.4 Neutrinos and the r-process\ 



Rather than describe the exotic effects of neutrino oscillations on the r-process, 
mentioned briefly above, we will examine standard-model effects that are never- 
theless quite interesting. The nuclear physics of this section neutrino-induced 
neutron spallation reactions - is also relevant to recently proposed supernova 
neutrino observatories such as OMNIS and LAND. In contrast to our first 
discussion of the ^-process in Sec. 4.2, it is apparent that neutrino effects 
could be much larger in the hot bubble r-process: the synthesis occurs much 
closer to the star than our Ne radius of 20,000 km: estimates are 600-1000 
km. The r-process is completed in about 10 seconds (when the temperature 
drops to about one billion degrees), but the neutrino flux is still significant as 
the r-process freezes out. The net result is that the "post-processing" neutrino 
fluence - the fluence that can alter the nuclear distribution after the r-process 
is completed - is about 100 times larger than that responsible for fluorine 
production in the Ne zone. Recalling that 1/300 of the nuclei in the Ne zone 
interacted with neutrinos, and remembering that the relevant neutrino- nucleus 
cross sections scale as A, one quickly sees that the probability of a r-process 
nucleus interacting with the neutrino flux is approximately unity. 

Because the hydrodynamic conditions of the r-process are highly uncer- 
tain, one way to attack this problem is to work backward in time. We know 
the final r-process distribution (what nature gives us) and we can calculate 
neutrino-nucleus interactions relatively well. Thus from the observed r-process 
distribution (including neutrino post-processing) we can work backward to find 
out what the r-process distribution looked like at the point of freeze-out. In 
Figs. 16 and 17, the "real" r-process distribution - that produced at freeze-out 
- is given by the dashed lines, while the solid lines show the effects of the neu- 
trino post-processing for a particular choice of fluence. The nuclear physics 
input into these calculations is precisely that previously described: GT and 
first-forbidden cross sections, with the responses centered at excitation ener- 
gies consistent with those found in ordinary, stable nuclei, taking into account 
the observed dependence on — Z\. 

One important aspect of the figures is that the mass shift is significant. 
This has to do with the fact that a 20 MeV excitation of a neutron-rich nucleus 
allows multiple neutrons ( ~ 5) to be emitted. (Remember we found that the 
binding energy of the last neutron in an r-process neutron- rich nuclei was about 
2-3 MeV under typical r-process conditions.) The second thing to notice is that 
the relative contribution of the neutrino process is particularly important in 
the "valleys" beneath the mass peaks: the reason is that the parents on the 
mass peak are abundant, and the valley daughters rare. In fact, it follows 
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Figure 16: Comparison of the r-proccss distribution that would result from the freeze-out 
abundances near the A ~ 130 mass peak (dashed line) to that where the effects of neutrino 
post-processing have been include (solid line). The fluence has been fixed by assuming that 
the A = 124-126 abundances arc entirely due to the ^-process. 

from this that the neutrino process effects can be dominant for precisely seven 
isotopes (Te, Re, etc.) lying in these valleys. Furthermore if an appropriate 
neutrino fluence is picked, these isotope abundances are produced perfectly 
(given the abundance errors). The fluences are 

N = 82 peak 0.031 • 10 51 ergs/(100km) 2 /flavor 
N = 126 peak 0.015 • 10 51 crgs/(100km) 2 /fiavor, 

values in fine agreement with those that would be found in a hot bubble r- 
process. So this is circumstantial but significant evidence that the material 
near the mass cut of a Type II supernova is the site of the r-process: there is 
a neutrino fingerprint. 

In the following sections we will briefly return to the effects of supernova 
density and electron fraction fluctuations and sterile neutrino mixings on the 
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Figure 17: As in Fig. 16, but for the A ~ 195 mass peak. The A = 183-187 abundances are 
entirely attributed to the ^-process. 



supernova r-process. 



6 MSW Exotica 

6.1 Spin-Flavor Oscillations 

In section 2.9 we briefly touched upon the effects of neutrino magnetic moment 
on neutrino propagation through matter. The chiral components of neutrino 
propagating in matter satisfy the equation!^ 
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where N e and N n are the electron and neutron densities in the media, respec- 
tively. As Eq. (124) demonstrates the chiral rotation is actually suppressed 
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in matter. It was soon realized E3, however, that combining the chiral preces- 
sion and the MSW flavor transformation by including non-diagonal magnetic 
moments circumvents this problem. One then should consider the evolution 
of left- and right-handed components of the electron and muon neutrinos to- 
gether: 
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For Dirac neutrinos the block matrices appearing in the evolution Hamiltonian 
in Eq. (125) are given by 
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and 



M = 



(128) 



In these equations B is the component of the magnetic field transverse to the 
direction of neutrino propagation. The last equation contains both diagonal 
and transition magnetic moments for neutrinos. For Majorana neutrinos di- 
agonal moments vanish and, since neutrinos are their own antiparticles, both 
left-handed and right-handed neutrinos must mix with the same mixing matrix. 
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Hence for Majorana neutrinos the evolution matrix takes the form 
^(2N e -N n ) *g 



r ^{2N e - N n ) ^-sva26 v n*B 



% sin 26 v % cos 29 v - ^N n —fi*B 



-liB -^(2N e -N n ) %■ 



sin 26,,, 



[iB ^sin20„ ^cos2^ + ^|iV n 

(129) 

These, evolution equations were numerically investigated in detail both for the 
Sun till and supernovaeEj. This was motivated by the Homestake data which 
suggested an anti-correlation between sunspot number and the capture rate 
of solar neutrinos. It is rather difficult to understand such a behavior in a 
theoretical framework since the most straightforward explanation for short- 
term time- variation of the neutrino flux is to assume the existence of a rather 
large neutrino magnetic moment, inconsistent with the bounds discussed in 
Section 3.3. Also Kamiokande and SuperKamiokande solar neutrino data do 
not show such time- variations. 

6.2 Neutrino Propagation in Stochastic Media 

In implementing the MSW solution to the solar neutrino problem one typically 
assumes that the electron density of the Sun is a monotonically decreasing funn- 
tion of the distance from the core and ignores potentially de-cohering effectsEj. 
To understand such effects one possibility is to study parametric changes in the 
density or the role of matter currents l21 In this regard, Loreti and Balantekin 
Eil considered neutrino propagation in stochastic media. They studied the sit- 
uation where the electron density in the medium has two components, one 
average component given by the Standard Solar Model or Supernova Model, 
etc. and one fluctuating component. Then the Hamiltonian in Eq. (53) takes 
the form 

H= \—^- cos 20+— G F (N e (r)+K(r))\a z + ^— S m29)a x . (130) 



GG 




Figure 18: Mean electron neutrino survival probability in the sun with fluctuations. The 
average electron density is given by the Standard Solar Model of Bahcall and Pinsonneault 

and sin 2 29 = 0.01. 



where one imposes for consistency 

(K(r)) = 0, (131) 

and a two-body correlation function 

(Ar(r)AT(r')) = 1 N e (r) N e (r') exp(-|r - r'\/r c ). (132) 

In the calculations of the Wisconsin group the fluctuations are typically taken 
to be subject to colored noise, i.e. higher order correlations 

f 12 ... = (N r e (n)N r e (r 2 )---) (133) 

are taken to be 

/l234 = /l2/34 + /l3./24 + /l4/23) (134) 

and so on. 

The mean survival probability Lj for the electron neutrino in the Sun is 
shown in Figure ^ where fluctuations are imposed on the average solar elec- 
tron density given by the Bahcall-Pinsonneault model. One notes that for 
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very large fluctuations complete flavor de-polarization should be achieved, i.e. 
the neutrino survival probability is 0.5, the same as the vacuum oscillation 
probability for long distances. To illustrate this behavior the results from the 
physically unrealistic case of 50% fluctuations are shown. Also the effect of the 
fluctuations is largest when the neutrino propagation in their absence is adia- 
batic. This scenario was applied to the neutrino convection in a core-collapse 
supernova where the adiabaticity_condition is satisfiedl£3. Similar results were 
also obtained by other authors Ea It may be possible to test solar matter 
density fluctuations at the BOREXINO detector currently under construction 
E3. Propagation of a neutrino with a_caagnetic moment in a random magnetic 
moment has also been investigated Also if the magnetic field in a po- 

larized medium has a domain structure with different strength and direction 



felt by the neutrinos 
Using the formalism 



in different domains, the modification of the potenti. 
due polarized electrons will have a random character E 
sketched above, it is possible to calculate not only the mean survival prob- 
ability, but also the variance, cr, of the—fluctuations to get a feeling for the 
distribution of the survival probabilities as illustrated in Figure [lj]. In these 
calculations the correlation length r is taken to be very small, of the order of 
10 km., to be consistent with the helioseismic observations of the sound speed 
(cf. Section 2.5). In the oppositeJimit of very large correlation lengths are 
very interesting result is obtained Ej, namely the averaged density matrix is 
given as an integral 



lim (p(r)) = -== / dxexp[-x 2 /(2^)]p(r,x), (135) 

reminiscent of the channel-coupling problem in nuclear physicsEl. Even though 
this limit is not appropriate to the solar fluctuations it may be applicable to a 
number of other astrophysical situations. 



6.3 Approximate Solutions and Neutrino Tomography 

Exact solutions for the neutrino propagation equations in matter exist only 
for a limited class of density profiles. For example, for those density profiles 
that satisfy an integrability condition called shapeJiivariance neutrino survival 
probability can be found with algebraic methods Ej. One can also use various 
approximation methods to solve Eq. (53\_ The standard semi-classical ap- 
proximation gives the adiabatic evolution 113. For a monotonimlly changing 
density profile supersymmetric uniform approximation yields Lil the hopping 
probability of Eq. (65) to be 

Phop = cxp(-7rfi), 
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Figure 19: Mean electron neutrino survival probability plus minus a in the sun with 
fluctuations. The average electron density is given by the Standard Solar Model of Bahcall 
and Pinsonneault and sin 2 29 = 0.01. Panels (a), (b), (c), and (d) correspond to an average 
fluctuation of 1%, 2%, 4%, and 8% respectively. 



n = ~ P dr [C 2 M ~ 2C(r) cos 29 v + l] 1/2 , (136) 

where Tq and r are the turning points (zeros) of the integrand. In this expres- 
sion we introduced the scaled density 

where N e is the number density of electrons in the medium. By analytic 
continuation, this complex integral is primarily sensitive to densities near the 
resonance point. The validity of this approximate expression is illustrated in 
Figure ^0|. As this figure illustrates the approximation breaks down in the 
extreme non-adiabatic limit (i.e., as Sm 2 — > 0). Hence it is referred to as the 
quasi-adiabatic approximation. 
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Figure 20: T he e lectron neutrino survival probability for the Sun. The solid line is calculated 
using Eq. (136). The dashed line is the exact (numerical) result. The dotted line is the 



linear Landau-Zener result. In the top figure, the lines are indistinguishable. An exponential 
density with parameters chosen to approximate the Sun was used. 



The near-exponential form of the density profile in the Sun motivates an 
expansion of the electron number density scale height, r s , in powers of density: 



(138) 



where prime denotes derivative with respect to r. In this expression a mi- 
nus sign is introduced because we assumed that density profile decreases a r 
increases. (For an exponential density profile, N e ~ e~ ax , only the n = 
term is present). To help assess the appropriateness of such an expansion the 
density scale height for the Sun calculated using the Standard Solar Model 
density profile is plotted in Figure ^l|. One observes that there is a significant 
deviation from a simple exponential profile over the entire Sun. However the 
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Figure 21: Electron number density scale height (cf. Eq. (138)) as a function of the radius for 
the Sun. The dashed line is the exponential fit over the whole Sun. The shaded are indicates 
where the small angle MSW resonance takes place for neutrinos with energies 5 < E < 15 

MeV. 



expansion of Eq. (138) needs to hold only in the MSW resonance region, in- 
dicated by the shaded area in the figure. Real-time counting detectors such 
as SuperKamiokande and Sudbury Neutrino Observatory, which can get infor- 
mation about energy spectra, are sensitive to neutrinos with energies greater 
than about 5 MeV (cf. Section 2.8). For the small angle solution (sin 29 ~ 0.01 
and 8m 2 = 5 x 10~ 6 eV 2 ), the resonance for a 5 MeV neutrino occurs at about 
0.35 R.0 and for a 15 MeV neutrino at about 0.45 (the shaded area in the 
figure) . In that region the density profile is approximately exponential and one 
expects tha t it should be sufficient to keep only a few terms in the expansion 
in Eq. (138) to represent the density profile of the Standard Solar Model. 

and using an integral 



Inserting the expansion of Eq. (138) into Eq. (|13(, 
representation of the Legendre functions, one obtains Ej 



n 



5m 
~2E 



|bo(l ~ cos20„) 
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Figure 22: (a) Spectrum distortion at SNO for the small-angle MSW solution {5m? ~ 
5 X 10" 6 eV 2 and sin 29 ~ 0.01). The solid line is the exact numerical solution. T he d ashed, 
dot-dashed, and dotted lines result from values of n up to 0, 1, and 2 in Eq. (139). The 
error bars on the exact numerical result correspond to two and five years of data collection. 
The dot-dot-dot-dashed line is the spectrum without MSW oscillations, normalized to the 
same total rate as with MSW oscillations. Note that on the scale of this figure the n = 1 
and n = 2 lines are not distinguishable from the exact ans wer, (b) The relative error arising 

from the use of Eq. (13U). 



+ 



E 



/ 5m 2 
{ \2V2G F E J 2n + 1 



[P w _i(cos20„) - P n+ i(cos20„)]} , (139) 



where P n is the Legendre polynomial of order n. The n — term in Eq. (139) 
represents the contribution of the exponential density profile alone. Eq. (139) 
directly connects an expansion of the logarithm of the hopping probability in 
powers of 1/2? to an expansion of the density scale height. That is, it provides a 
direct connection between N e (r) and P V {E V ). Eq. (139) provides a quick and 
accurate alternative to numerical integration of the MSW equation for any 
monotonically-changing density profile for a wide range of mixing parameters. 
The accuracy of the expansion of Eq. (139) is illustrated in Figure where 
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the spectrum distortion for the small angle MSW solution is plotted. One 
observes that for the Sun, where the density profile is nearly exponential in 
the MSW resonance region, the first two terms in the expansion provide an 
excellent approximation to the neutrino survival probability. 



6.4 Oscillations into Sterile Neutrinos 

In describing the MSW mechanism in Section 2.7 we assumed that the neutrino 
that mixes with the electron neutrino is an active neutrino (taken to be the 
muon neutrino for definiteness). We define an active neutrino to be one that 
couples to the Z, i.e. v e , Vfi, or v T . Data from the solar and atmospheric 
neutrino experiments, taken together with the intriguing data from the LSND 
experiment (cf. Section 8) may be viewed as an indication of a fourth neutrino 
flavor. Since the Z decay gives the number of active flavors to be three, if there 
is fourth flavor, it must be a sterile neutrino. (We define a sterile neutrino 
to be one that does not couple directly to Z and W, but may mix with those 
neutrinos that do). Many theories of neutrino mass predict the existence of 
sterile neutrinos. These are, however, generally rather heavy. It is not easy to 
build a model with sterile neutrinos light enough to mix with the active ones. 
A sterile neutrino, may have very interesting implications for the supernova 
r- process as wellE2l. 

MSW evolution equations describing the mixing of active and sterile neu- 
trinos take a different form than Eq. (53), In writing down Eq. (53) we used the 
fact that both active neutrinos acquire similar contributions to their effective 
masses from the neutral current interactions on neutrons. Such a term then 
only contributes to an overall phase which does not effect neutrino evolution. 
Since a sterile neutrino has neither charged- nor neutral-current interactions 
with ordinary hadrons we need to keep the neutron-term for the active flavor. 
Hence the MSW evolution equation takes the form 
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where 

<Pe(r) 



N-(r)~N+(r) 



N n (r) 



E - 8m 2 cos 26 v ) (141 



for the mixing of electron neutrinos (the plus sign on the right-hand side of 
the equation) or electron antineutrinos (the minus sign) with sterile neutrinos. 
Since matter-enhanced mixing of sterile and active neutrinos is typically used 
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for supernova environments in Eq. (141) we included the positron (N^(r)) as 



well as the electron(A f ~(r)) backgrounds. Note that in what follows, we take 
the sterile neutrino to be predominantly the heavier mass eigenstate. 

The mixing of muon and tau neutrinos with sterile neutrinos may be de- 
scribed similarly. The evolution Hamiltonian is as for the electron neutrino 
species, but with ip^ or ip T replacing ip e in Eq. ( |14C| ) as appropriate, where 

T (r) = -— (±V2 G F N n (r)E + 8m 2 cos 29 v ) ■ (142) 



As before, the + sign corresponds to neutrino mixing, and the — sign to an- 
tincutrino mixing. 

For astrophysical applications it is useful to express these equations in 
terms of the electron fraction. For a neutral medium we have Y p = Y e and 
Y n = 1 — y e , where Y p and Y n give the number of all protons or neutrons 
(free as well as those bound in nuclei), respectively, relative to baryons. The 
electron fraction Y e is given by 

N~(r) - N + (r) 

Y e (r) = e [ 1 , e [ ' . (143) 

y) Ne(r)-N+(r) + N n (r) V ' 



Inserting Eq. (143) into Eq. (141) one obtains the diagonal terms in the evo- 



, 3G F p(r) / 1\ 8m 2 „„ 
p e (r) = ± - if y ' Ye - r - -tft cos 26 v , (144) 



lution operator to be 
and 

^, T (r) = ±^r^ (Y e - 1) - % co S 20 v , (145) 

where p(r) is the matter density and tun is the nucleon mass. 

One should also point out that for the standard active-active MSW evolu- 
tion only either neutrinos or antineutrinos (depending the sign of 5m 2 ) undergo 



a resonance. The situation for active-sterile mixing is different. Eq. (144) in- 
dicates that, with appropriate neutrino parameters (i.e. 8m 2 > 0) and matter 
density, for Y e > 1/3 electron neutrinos and for Y e < 1/3 electron antineutri- 
nos can undergo an active-sterile MSW resonance. The possibility of matter- 
enhanced conversion of both iViL.and v^s can have interesting consequences 
in an astrophysical environment o. 

In the presence of neutrino fluxes ( "background" neutrino s) th e neutrino 
amplitude evolution Hamilton and the effective mass in Eq. ( |140| ) will have 
an additional term due to neutrino-neutrino neutral current forward exchange 
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scattering. In the case of active-active neutrino evolution, the neutrino back- 
ground,, because of flavor mixing, contributes to both diagonal and off-diagonal 
termsLij in the flavor basis amplitude evolution Hamiltonian, Eq. 
ever, for active-sterile mixing the off-diagonal terms are identically zerollfl. 

7 Neutrino Oscillation Constraints from the r-process 

As we discussed in Section 5 understanding neutrino transport in a supernova 
is an essential part of understanding supernova dynamics. Neutrino transport 
in a medium like supernova is a complicated process which needs to be treated 
numerically taking into account many different pieces of physics. In a core- 
collapse driven supernova, the inner core collapses subsonically, but the outer 
part of the core supersonically. At some point during the collapse, when the 
nuclear equation of state stiffens, the inner part of the core bounces, but the 
outer core continues falling in. The shock wave generated at the boundary 
loses its energy as it expands by dissociating material falling through it into 
free nucleons and alpha particles. For a large initial core mass, the shock wave 

fts stalled at ~ 200 to 500 km away from the center of the proto-neutron star 
. Meanwhile, the proto-neutron star, shrinking under its own gravity, loses 
energy by emitting neutrinos, which only interact weakly and can leak out on 
a relatively long diffusion time scale. One question to be investigated then 
is the possibility of regenerating the shock by neutrino heating. The density 
at the neutrinosphere is ~ I0 12 g cm -3 and the density at the position of the 
stalled shock isiza ~ 2 x 10 g cm . Writing the MSW resonance density in 
appropriate units: 

, / Sm 2 \ /MeVN /0.5\ , , x 

^ = i-3ixio 7 {-^){-e7)[tJ^ ' (146) 

one sees that, for small mixing angles, E u ~ 10 MeV, and cosmologically 
interesting 8m? ~ 1 — 10 4 eV 2 , there could be an MSW resonance point between 
the neutrinosphere and the stalled shock. 

Most neutrinos emitted from the core are produced by a neutral current 
process, and so the luminosities are approximately the same for all flavors. The 
energy spectra are approximately Fermi-Dirac with a zero chemical potential 
characterized by a neutrinosphere temperature. The v T , V T , i>u , F M interact with 
matter only via neutral current interactions. These decouple at relatively small 
radius and end up with somewhat high temperatures, about 8 MeV. The F e 's 
decouple at a larger radius because of the additional charged current interac- 
tions with the protons, and consequently have a somewhat lower temperature, 
about 5 MeV. Finally, since they undergo charged current interactions with 



140). r How- 
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more abundant neutrons, f e 's decouple at the largest radius and end up with 
the lowest temperature, about 3.5 to 4 MeV. An MSW resonance between the 
neutrinosphere and the stalled shock can then transform v T <-> i> e , cooling f T 's, 
but heating v^s. Since the interaction cross section of electron neutrinos with 
the matter in the stalled shock increases with—increasing energy, it may be 
possible to regenerate the shock. Fuller et al. l3 found that for small mixing 
angles between v T and v e one can get a 60% increase in the explosion energy. 

There is another implication of the v T and v e mixing in the supernovae. 
Supernovae are possible r-process sites (cf. Section 5.4), which requires a 
neutron-rich environment, i.e., the ratio of electrons to baryons, Y e , should he 
less than one half. Y e in the nucleosynthesis region is given approximately Ej 

by ^ | 

where A„ e „, etc. are the capture rates. Hence if > T Vcl then the medium 
is neutron-rich. As we discussed above, without matter-enhanced neutrino os- 
cillations, the neutrino temperatures satisfy the inequality T Vt > 2V e > T Vc . 
But the MSW effect, by heating v e and cooling v T can reverse the direction 
of inequality, making the medium proton-rich instead. Hence the existence of 
neutrino mass and mixings puts severe constraints on heavy-element nucleosyn- 
thesis in supernova. In turn requiring supernova to be an r-process site implies 
constraints on the neutrino parameters. These constraints are investigated in 
Ref. 54. 

If the supernova is an r-process site it is also desirable to have neutron to 
seed nucleus ratio > 100 in order that the heavier r-process species (i.e., those 
in the A = 195 peak) can be produced. This ratio is basically determined by 
three quantities: i) the expansion rate; ii) the electron, fraction Y e ; and iii) the 
entropy per baryon. Though different calculationsELlEj disagree on the value of 
the entropy in the neutrino-driven wind during the r-proccss nucleosynthesis, 
several models can produce values of these three parameters that yield a high 
enough neutron-to-seed nucleus ratio at freeze-out to effect a reasonable r- 
process. Unfortunately there are neutrino-induced processes operating during 
or immediately after freeze-out which can work to greatly reduce the neutron- 
to-seed nucleus ratio to the point where acceptable r-process nucleosynthesis in 
this site would be impossible. These neutrino-induced r-process destroyers are: 
i) neutrino neutral current spallation of alpha particles; and ii) the v e + n — ► 
p + e~ reaction accompanying the formation of alpha particles, also known as 
the "alpha effect." 

Meyer pointed out that previously neglected neutrino spallation reactions 
on the alpha particles tend to inhibit the r-process by allowing the assembly 
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of too many seed nucleill3. This process is especially effective at wrecking the 
r-process where the entropy is high. A simple steady-state wind model survey 
of the thermodynamics-parameters in neutrino-heated outnow was conducted 
by Qian and WoosleyllS. These authors concluded that the entropy in such 
models should be ~ 100/c per baryon, as opposed to Mayle and Wilson's model 
with an entropy of ~ 400fc per baryon. In turn, this result might argue against 
the effectiveness of neutrino-induced alpha particle spallation in lowering the 
neutron to seed nucleus ratio. However, lower entropies in general imply a lower 
value of this ratio since there will be more seed nuclei in these conditions. At 
best, the neutron to seed nucleus ratios obtained in lower entropy models are 
marginal for the production of the neutron-rich r-process species l 79 lrH. 

The alpha effect occurs at the epoch of alpha particle formation. As the 
temperature drops, essentially all the protons and most of the neutrons in 
the ejecta lock themselves into alpha particles which have a large binding 
energy. This phenomenon ultimately will tend to push the electron fraction 
higher, towards Y e = 0.5. The increase in Y e comes about because protons 
produced by electron neutrino capture on neutrons will in turn capture more 
neutrons to bind into alpha, particles, reducing the number of free neutrons 
available for the r-process c3. This effect has been shosscn to be the biggest 
impediment to achieving an acceptable r-process yield E3. Matter-enhanced 
neutrino transformation between electron neutrinos and other active species 
worsen this problem as it tends to increase electron neutrino energies. Hence 
one would like to avoid active-active neutrino mixing. Indeed this is how the 
bounds of Ref. 54 on neutrino parameters were obtained. 

One way to avoid or reduce the efficacy of the alpha effect is to reduce the 
flux of electron neutrinos at some point above the surface of the neutron star. 
However, in models of the neutrino-driven wind a large flux of electron neutri- 
nos is required to lift the material off the surface of the neutron star. In fact 
since nucleons are gravitationally bound by about ~ 100 MeV near the surface 
of the neutron star, and since each neutrino has an energy ~ 10 MeV, each nu- 
cleon must suffer some ~ 10 neutrino interactions to be ejected to infinity. So 
if we are to reduce the v e flux we must do so only at relatively large radius, so 
that effective neutrino heating already can have occurred. This can be achieved 
by the matter-enhanced active-sterile neutrino transformation in the v e v s 
and v e # v s channels. In such a schemeEj the lightest sterile neutrino would 
be heavier than the v e and split from it by a vacuum mass-squared difference of 
3 e V 2 < Sm^s < 70 eV 2 with vacuum mixing angle sin 2 29 es > 10~ 4 . Further 
details for this scenario is given in Ref. 52. A similar mechanism exploiting 
different active-sterile mixing channels is presented in Ref. 83. 
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8 Other Neutrino Oscillation Experiments 



An introduction to neutrino oscillations in vacuum was given in Section 2.6. In 
this section we briefly summarize some recent results from several neutrino os- 
cillation experiments. From Eq. (47) one can write the appearance probability 
of the other flavor 

P = sin 2 26 v sin 2 (1.275m 2 L/E). (148) 

In Eq. (148) Sm 2 is measured in eV 2 and L/E in m/MeV. Neutrino oscillation 
experiments are somewhat arbitrarily divided into two classes: short-baseline 
and long-baseline .]] As Eq. ( 148 ) indicates, the longer is the baseline, L, the 
more sensitive the experiment is to smaller values of Sm 2 . 

The Liquid Scintillator Neutrino Detector (LSND) collaboration at Los 
Alamos National Laboratory reported a significant "oscillation-like" excessEa. 
This detector sits 30 meters away from the beam dump at the old Los Alamos 
Meson Physics Facility (LAMPF). The experimental apparatus is designed to 
produce a beam of muon antincutrinos with as little contamination as possi- 
ble from the electron antineutrinos. If these muon antincutrinos oscillate into 
electron antineutrinos, such formed electron antineutrinos would interact with 
protons in the detector, creating a positron and a neutron. This neutron, after 
some time, binds with a proton to form a deuteron, giving a photon with a 
characteristic energy of 2.2 MeV. The experiment observes these photons as 
well as the positron's Cerenkov track. When they identify both signatures 
together, the LSND group first reported seeing nine events versus an expected 
background of two events coming from the electron antineutrinos from sources 
other than the muon antineutrino oscillation. This signal received a lot of 
attention and a fair amount of criticism. Indeed, a dissenting member of the 
LSND team has performed a data analysi^of his own which finds no posi- 
tive signal above the expected backgroundEJ. The LSND result is unlikely J-d 
be a statistical fluctuation. The excess persists after analyzing more dataEj, 
including /x decay in flight, which has different systematics and backgrounds 
from the decay at rest analysis. However, the KARMEN collaboration, carry- 
ing out a similar (but not identical) experiment at Rutherford Laboratory in 
England, reported no evidence for neutrinn. oscillations in a parameter space 
which largely overlaps with that of LSND EZL 

Recently CHOOZ, a long-baseline reactor-neutrino vacuum-oscillation ex- 
periment, announced E3 first results. This experiment, operating in a reactor 
in France detects electron antincutrinos by a liquid scintillation calorimeter lo- 
cated at a distance of about 1 km from the source. From the statistical agree- 
ment between detected and expected neutrino event rates, they find (at 90in 

c Atmospheric neutrino experiments can be considered very-long baseline experiments. 
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the electron antincutrino disappearance mode for the parameter region given 
approximately by 8m 2 > 0.9- 10 -3 eV 2 for maximum mixing and sin 2 29 > 0.18 
for large values of 5m 2 . The measurement of the atmospheric electron and 
muon neutrino zenith angle distributions at SuperKamiokande, taken together 
with the CHOOZ data presents a strong evidence for the mixing of muon 
neutrinos with either tau neutrinos or sterile neutrinos. 

This work was supported in part by the National Science Foundation, the 
US Department of Energy, and the University of Wisconsin Research Commit- 
tee with funds granted by the Wisconsin Alumni Research Foundation. 
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